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FIELD THEORY FOR TRAPPED ATOMIC GASES 



H.T.C. Stoof 



Abstract 

In this course we give a selfcontained introduction to the quan- 
tum field theory for trapped atomic gases, using functional meth- 
ods throughout. We consider both equilibrium and nonequilibrium 
phenomena. In the equilibrium case, we first derive the appropriate 
Hartree-Fock theory for the properties of the gas in the normal phase. 
We then turn our attention to the properties of the gas in the super- 
fluid phase, and present a microscopic derivation of the Bogoliubov 
and Popov theories of Bose-Einstein condensation and the Bardeen- 
Cooper-Schriefler theory of superconductivity. The former are ap- 
plicable to trapped bosonic gases such as rubidium, lithium, sodium 
and hydrogen, and the latter in particular to the fermionic isotope of 
atomic lithium. In the nonequilibrium case, we discuss various top- 
ics for which a field-theoretical approach is especially suited, because 
they involve physics that is not contained in the Gross-Pitaevskii 
equation. Examples are quantum kinetic theory, the growth and 
collapse of a Bose condensate, the phase dynamics of bosonic and 
fermionic superfluids, and the coUisionless collective modes of a Bose 
gas below the critical temperature. 

1 Introduction 

An important trend in the condensed matter physics of the last two decades, 
has been the use of advanced field-theoretical methods to discuss various 
subtle and fundamental properties of interacting many-particle systems at 
low temperatures. There are several reasons for this trend. The first reason 
is of course, that a traditional topic in statistical and condensed matter 
physics is the study of phase transitions and critical phenomena, for which 
the universal properties are independent of the microscopic details of the 
system and can therefore be determined by a field theory describing only the 
large-scale properties of the system of interest. Since the latter are usually 
solely determined by symmetry considerations, this has led to the important 
concept of spontaneous symmetry breaking, which has turned out to be not 
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only highly successful in condensed matter physics, but also in high-energy 
physics and in particular in the Standard Model of elementary particles [Q . 

A second reason is that soon after the development of the renormaliza- 
tion group methods for critical phenomena l2l , it was realized that the same 
methods can in fact be used to describe the large-scale properties of many- 
particle systems at any temperature and not only near the critical one. 
Moreover, application of the renormalization group ideas does not only lead 
to an understanding of the static behaviour but also of the dynamical prop- 
erties H. As a result quantum field theories can be used to conveniently 
determine the dynamics of many-particle systems close to equilibrium, i.e., 
for example the linear hydrodynamical equations of motion. In addition, it 
can even be used in highly nonequilibrium situations where in general also 
nonlinearities play an important role. This feature, that also the dynamics 
of the system can be captured by a quantum field theory, has for instance 
led in recent years to the study of so-called quantum phase transitions [H . 

Finally, the importance of field-theoretical methods in condensed matter 
physics is associated with the observation that also the effects of imperfec- 
tions, i.e., disorder, can be treated in this way B. Apart from the tech- 
nological importance of disorder, for example for superconducting magnets, 
disorder leads also to fundamentally new physics such as the phenomenon 
of localization [pj. In quantum Hall systems, a combination of disorder and 
interaction effects give rise to the realization of various peculiar quantum 
fluids with fractionally charged excitations M . The application of field theo- 
ries has proven to be highly successful in this case and has led to a theory of 
the quantum Hall effect in terms of edge states that form a chiral Luttinger 
liquid 1^. In mesoscopic physics, the study of disorder in small electronic 
structures has resulted in the so-called random matrix theory ||9[ , which has 
also been of much use in the study of the quantization of classically chaotic 
systems. 

In this course we aim to show that quantum field theory is also very 
convenient for obtaining a detailed understanding of the equilibrium and 
nonequilibrium properties of trapped atomic gases. After the first observa- 
tions of Bose-Einstein condensation in 1995 ||l(]|-[l^, degenerate atomic gases 
have received again a great deal of attention and are presently the main 
subject of study of a large number of experimental and theoretical groups 
around the world. The reason for all this excitement is, first of all, that 
Bose-Einstein condensation has never before been observed experimentally 
in a clear-cut manner, even though this phenomenon was already predicted 
by Einstein in 1925 [|l3[. Second, it is of fundamental interest because it is 
the only phase transition that occurs also in the absense of interactions and 
is, therefore, the textbook example for the use of statistical-physics meth- 
ods. Finally, the goal of achieving Bose-Einstein condensation turned out 
to be much more difficult than anticipated at first. Just before the break- 



H.T.C. Stoof: Field theory for trapped atomic gases 



through in 1995, it had even acquired the nature of a 'quest for the holy 
grail', since the pioneering experiments were already performed in 1980 [ [L4[ . 
From a theoretical point of view, a quantitative understanding of the 
experiments with cold atomic gases requires that we take into account the 
following two effects. First, the gas is magnetically trapped in an, usually 
axially symmetric, harmonic oscillator potential. This is necessary because, 
in order to obtain the required record low temperatures of 1 — 100 nK, the 
gas cannot be allowed to have any contact with material walls. Second, 
the atoms of the gas interact with each other, which in general dramatically 
affects the behaviour of the condensed gas, even at the relevant low densities 
of 10^^ — lO^'' cm~'^. An accurate description of these degenerate gases 
thus requires the solution of a highly inhomogeneous many-body problem. 
The theoretical challenge posed by these new quantum systems lies therein 
that the density of the gas is sufficiently small that it should be possible 
to accurately solve this many-body problem from first principles and to 
compare the outcome of the theory directly with experiment, i.e., without 
any adjustable parameters. In our opinion, quantum field theory is the most 
simple way in which we are able to meet this challenge. 

2 Equilibrium field theory 

We start our development of the quantum field theory of trapped atomic 
gases by considering first the equilibrium properties of these gases. We 
consider both Bose and Fermi gases, and the ultimate aim of this section 
is to arrive, for both cases, at an accurate description of the normal and 
superfluid phases of the gas. Although mixtures of Bose and Fermi gases are 
also of current interest, we do not consider them explicitly here, because they 
can be treated by a straightforward generalization of the theory. In section 
we then turn to the nonequilibrium properties, which are perhaps the most 
interesting and certainly the least understood at present. The reason for 
organizing the course in this way, is that the development of the equilibrium 
theory gives us an opportunity to introduce all the necessary tools that 
are required for a treatment of the more complicated nonequilibrium case. 
In particular, we present in detail the way in which we can make use of 
functional methods. To make also a connection with the more familiar 
operator language, however, we first briefly summarize the outcome of the 
second quantization formalism. 

2.1 Second quantization 

The atoms of interest to us have internal degrees of freedom due to the 
electron and nuclear spins. In principle this is very important, because it 
gives the atom a magnetic moment, which is used to trap the atoms in a 
magnetic field minumum. During this course, however, we restrict ourselves 
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to atomic gases that are a mixture of at most two hyperfine states. Without 
loss of generality, we can then suppose to have N identical atoms with mass 
m and effective spin s in an external potential V°^{x). As a result, the 
time-dependent Schrodinger equation we have to solve is 

zn||*(t))=iJ|vI/(0), (2.1) 

where the hamiltonian is 

^ = E 1^ + ^'''(^0 - 7S. • B^M + - ^ y(x, - X,) , (2.2) 

[Xj,pj]_ = ih, and all other commutators of the positions and momenta 
vanish. The first term in the right-hand side is the sum of the one-particle 
hamiltonians, which includes an effective Zeeman interaction that accounts 
for a possible difference in the hyperfine energies. The second term repre- 
sents the interactions. For simplicity, we have assumed that the interaction 
Vijci — Xj) is independent of the hyperfine states of the atoms i and j. 
This is in general not justified for realistic atomic gases, but is valid for the 



specific applications that we have in mind. Moreover, in section 2.6 we also 
discuss the general case. Finally, we have also neglected possible three-body 
forces. This is a result of the fact that we are interested in dilute quantum 
gases, for which it is highly improbable for three atoms to simultaneously 
interact with each other. 

Without interactions the eigenstates are, of course, given by the states 
{|ni, q;i)i (g) |n2, 0:2)2 ... \nN,aN)N}, where the specific quantum state 
for each atom is exactly known. Here n = (rix, Uy, n^) and the nonnegative 
integers n^,, ny, and Uz denote the three quantum numbers that are re- 
quired to specify the one-particle eigenstates in the external potential. The 
wave functions and energies of these eigenstates are Xn(x) = (x|n) and e„, 
respectively, and are found from the time-independent Schrodinger equation 



2m 



+ V-^^) - e„ xn(x) = . (2.3) 



In addition, the internal state \a) is a shorthand notation for \s,ms). The 
many-body wave function |$(t)), however, has to be symmetric or anti- 
symmetric under permutations for bosonic or fermionic atoms, respectively. 
Therefore it is more convenient to use a properly (anti)symmetrized version 
of the above basis, i.e., the states |{-/Vn,Q}) with the occupation numbers 
^^11,0 = 0, 1, 2, . . . , 00 for bosons and N^^a = 0, 1 for fcrmions. The Hilbcrt 
space of all these states, without the constraint N = J2n a -^n.a, is known 
as the Fock space. 
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Clearly, the many-body wave function |^(i)) can be expanded in this 
basis as 

l*W>= E *({^n,4,0l{A^n,4), (2.4) 

where ^E'({-/Vn,a}, t) is the amplitude for the gas to be in state |{A^n,a}) at 
time t. In this basis the Schrodinger equation becomes 

^h^^^{{Nn,^},t)^ J2 {{Nn,c.}\H\{KJ)^{{KJ,t) , (2.5) 

which shows that we need the matrix elements of the hamiltonian between 
different states in the Fock space. To calculate these most easily we intro- 
duce so-called annihilation operators V'n.a by 



^n,a\...,Nn,a,...) = V^n^a | ■ ■ • , iV„,„ - 1, . . .) , (2.6) 

from which it follows that the creation operators ^pl^ ^ obey 

i'ij ■■■: ^n,a, ■ • ■) = ^1 ± Nn,c.\ . . . , A^„,a + 1, . . .) . (2.7) 

As a result, we see that the operator -0^ aV'n.a counts the number of atoms 
in the state |n,Q;), i.e., 

IpUnM . . . , A^„,a, . . .) = Nn,a\ ■ ■ • , A^„,a, ■ ■ •) ■ (2.8) 

We have also that [ipn,a,''Pn',a']^ — ['4'n a^^li' a'^T ~ ^ ^^"-^ most impor- 
tantly that 

[V'n,a,'0i/ o'It = Sn,n'Sa,a' ■ (2.9) 



In equations (2.7) and (2.9) the upper sign refers to bosons and the lower 
to fermions. This will be true throughout the course. 

From these results we can now easily show, first of all, that the basis in 
the Fock space is given by 

|{iVn,J)=n 7^ , |0)^ (2-10) 

with |0) the vacuum state containing no atoms. Second, the hamiltonian is 

H = Y,^n^c.i^lJn.c. (2.11) 

+ 2 ^ ^ Ki,n';m,m''0i,aV'I',a'V'm',a'V'm,a , 
Oi.a' n,n' m.m' 
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a) n,a 



n',a' b) 





in,a 



m ,a 



Fig. 1. Diagrammatic representation of tlie interaction terms in a) equation (2.11) 
and b) equation (2.13). 



with en,Q = Cn — jnisB — en + ea the one-particle energies. Furthermore, 
Jd^Jd^' x;(x)x;.(x')T^(x-x')Xn.(x)xm'(x') (2.12) 

is the amplitude for a collision which scatters two atoms out of the states 
|m, a) and |m', a') into the states |n, a) and |n', a'), as schematically shown 
in figure 1(a). 

Introducing the field operators V'a(x) = X^n V'n.QXn(x) and V'i(x) = 
Sn ''/'n aXnlx), that annihilate and create atoms in the spin state |a) at 
position X respectively, we can rewrite this result into 

H = ^ /rixV-Kx) 

+ ^E/'^^/'^^' ^i(x)V^I'(x')nx-x')^«'(x')^„(x) . 

a, a' 

Note that, due to the (anti) commutation relations of the creation and an- 
nihilation operators and the orthogonality of the wave functions x„(x), 
the field operators obey [^c«(x), ^q,'(x')]^ — [■)/Jt (x), ^]j, (x')]ip = and 
[■0a(x),-0]^, (x')];p = S(x — x')Sa^a'- Moreovcr, it is also important for the 
following to note that the number operator is 



2m 



+ y«'^(x) + e„U„(x) (2.13) 



^ = I] ^i,a V-n.c = 51 / ^^ ^" (^)^" (^) 



(2.14) 



and, similarly, that the effective total spin operator is 



S= $] V-LialslaO^n.a' = ^ /"dx V-t(x)(a|s|a')V'a'(x) . (2.15) 

n rv rv' rv rv' 
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The density of atoms in the state \a) is thus simply ria(x) — ip'^{x.)ipa{'x.). 
This completes our brief discussion of the second quantization formal- 
ism. In principle, we could now proceed to develop an operator formulation 
of the quantum field theory of interest to us. Since the experimentally 
most important observables can be expressed as appropriate products of 
the field operators, as we have just seen, this would essentially amount 
to the study of the (imaginary) time evolution of the Heisenberg opera- 
tor V'(x,r) = e(^-''*)^/''-0c,(x)e-(^-''^)^/'' at a fixed chemical potential 
fi [^. Put differently, the desired quantum field theory would be defined 
by the Heisenberg equation of motion 

?i— V-„(x,t) = [ff-M7V,v3a(x,r)]_ , (2.16) 

and we would need to solve this equation in a sufficiently accurate ap- 
proximation. As mentioned previously, however, we here want to develop 
Feynman's 'path-integral' formulation of the problem, which will turn out 
to be much more convenient for our purposes. To do so in a manner that 
is the same for both bosonic and fcrnfionic atonfic gases, we first need to 
introduce some mathematical background. 



2.2 Grassmann variables and coherent states 

We have seen that in the case of fermions, we need to make use of anti- 
commuting creation and annihilation operators. This automatically builds 
in the Pauh principle in the theory, because it implies that {ipl^ a)^|0) = 
and thus that the occupation numbers Nn,a are restricted to be either or 
1. For reasons that will become clear in a moment, it is in that case also 
convenient to introduce anticommuting complex numbers or Grassmann 
variables. The simplest example is to have two such Grassmann variables, 
say (f) and </)*. The set {1, (j), (p* , (/'*0}, and linear combinations thereof with 
complex coefficients, form then a so-called Grassmann algebra. 

By definition we have [(/>, 0]+ — [0, (/)*]+ — [0*, (/)*]+ = and thus in 
particular cj)'^ = 0*^ ~ 0. Therefore, the above set is indeed complete. The 
complex conjugation in this algebra is defined by {(/))* = cj)*, {4>*)* = 4>, 
and {(t>*(t>)* = {4>)*{4>*)* — (()*(t>- Moreover, we can also define an analytic 
function on this algebra by 

A(0*, 0) = an + ai2(t> + aai^* + a22(j)* (j) ■ (2.17) 

As a result, it is natural to define also a differentiation by 

^A((/)*,<^) = ai2-a22<^* . (2.18) 

ocp 
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To be m.ore precise this is in fact a left differentiation and the minus sign 
occurs, because we need to permute (f>* and before we can differentiate 
with respect to (j). So, similarly, we have 

—A{^\<t>)=a2i+a22(t> (2.19) 

and 

^^(^%^) = -^^(^%^) = --- (2.20) 

Next we also need integrations over these Grassmann variables. Note 
that since (/)^ = we have only two possible integrals, namely J d(j) 1 and 
J d(j) (j). We define these by 



^1 = (2.21) 

and 

d(p (p^l . (2.22) 

This means that integration is equivalent to differentation. The main reason 
for the above definition is that we want the integration to obey the usual 
rules of partial integration. In particular, this implies that 

aF(0) 



= , (2.23) 

for any function F((j)) = /i + /20. It is clear that this condition requires that 
J d(j) 1 = 0. The result oi J d(j) (j) is then solely a question of normalization. 
It turns out that we are primarily interested in integrals of the form 



d(l3*d<j)A{(t)*,(j})= / d0*d(/i(aii + ai29i+a2i(/)*+a22 0*0) = -a22 (2.24) 



as we will see in section 2.3 



Clearly, we can now also consider the Grassmann algebra based on the 
variables 0„ and (/)* with n = l,2,...,oo. What we will need in the following 
are gaussian integrals over these variables. It is not difficult to show with 
the above definitions that 

/ n d^ndc^n exp I - ^ C^n,«"^«' > = dct^ = e^["^^l . (2.25) 

n I n,n' J 

If the variables (j)n and (/)* were just ordinary complex numbers we would 
in contrast have the result 



n I n,n' I 



= e-^rlinA] (2.26) 

det^ 
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These last two results will be used many times in the following. 

One immediate use of these Grassmann variables is that we can now con- 
sider eigenstates of the annihilation operator 4'n,a, also when we are dealing 
with fermions. These eigenstates are called coherent states. Consider the 
state 

!</)„,„) = (1 - (/.„,„< J|0) = cxp {-(/.„,« V^i,„} |0) , (2.27) 

where (/)n,Q is a Grassmann variable that also anticommutcs with the cre- 
ation and annihilation operators in our Fock space. Clearly we have that 

V'n,a|'/'„.c) = l/'„,a|0) + (/)„,a V'n.aV'i,^ |0) (2.28) 

= 4>ii,a\0) = 0n,a(l - 4'ii,a'4'n,a)\^) = 4>ii,a\4'ii,a) , 

SO \(t>ii,a) is indeed an eigenstate of V'n,Q with eigenvalue 0n,a- In general 
we can now make the states 

|(/))=cxp|-^(/)„,„^t^J|0) (2.29) 

that obey 'il'n,a\4') = 4'n.a\(t')- Introducing the Grassmann-valued field 
^a(x) = X]ii'^n.aXn(x), the latter two relations can be rewritten as 

|(/))=exp|-5]|dx</)„(x)^l(x)l|0) (2.30) 

and tpai^^M) = <t>a{x)\4>)- 

It is important to note that these coherent states are not orthonormal. 
In contrast, we find that 

i^m - 11(01(1 - ^„,aC, j(i - <i^LJijm (2.31) 

= 11(1 + Ca'^n.a) = exp J ^ K,a^n,a [ 
n.a \ n,a ) 

= exp I ^ I dx C (x)(/)^ (x) I ^ e(-^l-^') . 
Nevertheless they obey a closure relation, as can be seen explicitly from 
fYl'^<Plo.d<f>n.,'.e^pl-J2<f>lc.Ko.\\^){<f>\ (2.32) 

n.a l^ n,a ) 

n,Q "^ 
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We write this from now on simply as 

:-('^l«|0)(,/,| = i . (2.33) 



The interesting observation at this point is that essentially the same 
formulaes also hold for bosons ||lq,0. We have only a minus sign difference 
in 

|(/))=exp|^y'dx0„(x)V;t(x)l|o) , (2.34) 

but then it can still be easily shown that ■0q(x)|0) = 0q,(x)|(/i). A convenient 
way to do so is, for instance, by noting that the commutation relation 
[V;q(x),-0^, (x')]_ = (5(x — x.')5a,a' imphes that V'q(x) acts as (5/5?/;^ (x) 
on these states. Furthermore, the same expressions hold for the overlap 
{4>\4'') and the closure relation, if we define the integration measure by 
J d[(f)*]d[(l)] = /nn,a'^'/'n,a#n,a/(27ri) in this case. 

Summarizing, we have thus found for bosons and fermions that 

|0)=exp|±^y'dx0„(x)V't(x)l|O) , (2.35) 

= e(^l^') , (2.36) 



and 

""^|0)((/.| = i. (2.37) 

The last ingredient we need is that in terms of these coherent states the 
trace of an operator O over the Fock space can be expressed as 

Tr[0] = I d[<j)*]d[<j)\ e-f-^l-^' (±01010) . (2.38) 

The minus sign in the fermionic case is easily understood from the fact that 
we then have ({iV„,a}|0)(</'|{A^n,a}) = (-0|{iV„,a})({A^„,a}|0), due to the 
anticommuting nature of the Grassmann variables. After this mathematical 
interlude we can now return to physics and to our goal of arriving at a unified 
treatment of bosonic and fermionic atomic gases. 

2.3 Functional integrals 

We are interested in determining the equilibrium properties of a trapped 
atomic gas at some temperature T . From statistical physics we know that all 
these properties can be obtained from the grand-canonical partition function 



Tr 



g-/3(H-MA') 



(2.39) 
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"M 



--np 



M-l 



T0=0 



Fig. 2. Illustration of the slicing of the imaginary time interval [0, hp] needed for 
Feynman's path-integral formulation of the partition function. 



where (3 = l/fcsT and /i the chemical potential. We thus need to evaluate 
this quantity. As mentioned previously, we want to do so by making use of 
quantum field theory and, in addition, of Feynman's path-integral approach 
to quantum mechanics. We thus want to write the partition function as a 
functional integral over time-dependent fields 0ct(x, r), just like the partition 
function of a single particle in an external potential can be written as a 
functional integral over all time-dependent paths x(t). This can indeed be 
achieved with the coherent states that we have introduced in the previous 
section. 



We start with using our formula (2.38) for the trace of an operator. 



Z= / d[^*]d[cj)] e-("*l"*)(±0|e-'5(^-^^)| 



(2.40) 



and observe that we are thus faced with the task of calculating the matrix 
elements ((/>M|e-''(^~''^)|?io) with (/)o(x) = <?!)(x) and 0;^(x) = ±(/)*(x). 
This is difficult in general but can be simplified in the following way. We first 
realize that the operator e~^'^^~^^^ is identical to the quantum mechanical 
evolution operator U{t,Q) = (,-KH-tJ.N)t/h evaluated at i = -ihp. Put 
differently, we thus want to calculate the matrix elements of the imaginary- 
time evolution operator U{—iT, 0) for r = hp. To do so, we next split the 
time interval [0, h(3] into M pieces, with r,„ — nifi(3/M and m = 0, 1, . . . , M. 
So At = hf3/M . The procedure is summarized in figure 2. 

At each intermediate time Tm we then apply a closure relation of the 
coherent states. This gives 






(2.41) 



^-AT{H-nN)/h 



\4>n 



Now we can use that in the limit M ^ oo we only need to know the latter 
matrix elements up to order At, because terms of order (Ar)^ lead to 
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corrections of order M(At)^ cx 1/M, which vanish in that hmit. Hence 
(</'™|e-^"(^-^^^/''|0™-i) ~ (0™|1 - ^t{H - ^lN)|h\ct^.m-l) (2.42) 

with the grand-canonical hamiltonian functional resulting from equations 



( |2.13[ ) and ( |2.14D equal to 

2t-t2 



i^r , 0] - E/ '^^ '^"(^) {-^ + ^"W + eo " /^l 0o(x) (2.43) 

+ ^E/«'^/«'^' C(x)C'(x')^(x-x')<^c.'(x')0a(x) , 
a, a' 

since Va(x)|0) = </)« (x) | </)) and ('?!)|?Al(x) = (<?!)|0;(x). 

Thus, neglecting terms of order (Ar)^, we can again exponentiate the 



right-hand side of equation (2.42), which leads in first instance to 



and therefore for the desired matrix clement of the imaginary-time evolution 
operator to 

{<i>M\e-f"^"~^^^M = / ( n dW.n]d[<l>m] e-(*"l^™) j (2.45) 



This can then be manipulated into the suggestive form 

(</>M|e-'^(^-^^)|^o) = e^^-'^-) J ( n dirraWm] j (2.46) 

J 1 Y^ A fi-i't>ni\(l>7n)- {(l>7n\(l>7n-l) , „r,* , 

xexp<^ -- 2^ At Ih — + H[(f)^,<j>rn-i_ 

Taking now the continuum limit Af — > oo and putting 0„j = (J){t„i), we 
find that 

(^M|e-''(*-''*)|0o) (2.47) 
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where the so-called Euclidean action is 




(2.48) 
*(x,r)?i|-</.„(x,r)+iJ[<^*(r),</.(r)]l . 

This is essentially the desired functional integral over the complex fields 
^a(x, r) with the boundary conditions 0q(x, 0) = (/)o;a(x) and (^*(x, ?i/3) = 
0^.Q,(x). It is precisely the field theory analogue of the Feynman path 
integral. To obtain the partition function we only need to put (/'o(x) equal 
to ±(/)jvf (x) and perform a last integration over (/)m(x) and (/)^^(x). It then 
finally becomes 

Z = f d[<j)*]d[(t)]e~^^'f''''f'^^'' (2.49) 



with the boundary conditions 0c((x, ?i/3) = ±(/)q,(x, 0), i.e., the fields are 
periodic in [0, hp] for bosons and antiperiodic for fermions. Note that in 



equation (2.47) we have used the same notation for the integration measure 



as in equation (2.49), although there is in principle one more integration in 
the expression for the partition function. The difference is in the continuum 
limit accounted for in the boundary conditions, which are in practice usually 
left implicite. Having arrived at an exact identity between the partition 
function and a functional integral, we are now going to familiarize ourselves 
with this identity, and with how to perform functional integrals in general, 
by considering the ideal quantum gases. 

2.4 Ideal quantum gases 

Since the partition functions Zq of the ideal quantum gases are known ex- 
actly, they are ideal test cases for our field-theoretical methods. Moreover, 
a thorough knowledge of the ideal quantum gases is also an important first 
step in understanding experiments with trapped atomic gases, because the 
effects of the interatomic interaction can essentially be included perturba- 
tively. How this perturbation theory is performed is discussed in detail in 



section 2.5, but before we can do that we need to understand the noninter- 



acting gases first. In that case, we have 

S'o [</)*,(/.] =^ / dT dx (2.50) 





2v72 



dr 2m 
and the partition function is a gaussian integral, which explains why we 



X 0;(x,T)^n— -__ + V^-(x)+6„-M^0„(x,t) 



were interested in gaussian integrals in section 2.2. It will be illustrative to 



evaluate this partition function in three different ways. 
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2.4.1 Scmiclassical method 

In the first way, we perform the evaluation of the trace involved in the 
definition of the partition function at the end of the calculation, and start by 
considering the matrix element {±(j)\e~^^'^~^^' \(j)) as the functional integral 

'^(O)=0 Vn.a / In.a J 



{1 f^'^ ( d 



Xexp{-- dT}6*{T)[h—+en,a-^^)4'nA'^)} ' 



which is the product for each n and a of the path integral 

d[0*]dM exp I 0* {h(3)cl>{hf3) - i y dTcf,*{T)(n^+e- fi) </)(t) I 

with the boundary conditions 0(0) = 4> and 4>{h(]) = zL(j>* . It clearly corre- 
sponds to the matrix element (±0|e~''''^~^)'^ ^\(j)) for one particular value 
of n and a. We calculate this matrix element by performing a shift in the 
integration variables, i.e., (/)(t) = 0ci(t) + ^(r) and (f)*{T) ~ (pliiT) + S.*{t), 
where (pciiT) obeys the 'classical' equations of motion 



sso[r,<i>] 



6(1)* (t) 



?i^+e-A*)0ci(T) = O (2.51) 



and similarly for (/)*i(t). The solutions with the correct boundary solutions 
are <?!)ci(t) = (j)e-^'-'''>''/^ and 0*i(t) = ±(j,*e^'-i''>^^-''l^'>/'', leading to the 
path integral 

exp{±e-^('-^V*0} 

X J d[Cm] exp I - i £^ dT C {r)(n^+e- ^^ Hr) \ 

with the boundary conditions ^*{'h(3) = ^(0) = 0. This means that the last 
path integral is just equal to {0\e~^^'^~'^'>'^ ^\0) = 1 and that the desired 
result is just the prefactor. 

Substituting this for each value of n and a, we apparently have 

(±0|e-^(^-^*) !</)) = exp |±^e-''(^"--^)C,a0n,4 (2.52) 
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and the partition function becomes the product 

= n(lTe-'^(^--''))Ti=exp|TEMlTe-^(^--''))| . 

n.a l^ n.a ) 

This is the correct resuh, because from the usual thermodynamic identity 
{N) = dlnZo/diP^i) we find 

W - E .,<', ^ , ' (2-54) 



e,3(e„,.-p) ^ 1 ' 
n,a 



as desired. 



2.4.2 Matsubara expansion 

The second way is easier and more common in practice. We immediately 
start with 



Za= d[(t)*]d[(j)] e-^ol'^'.'^l/'i (2.55) 

and incorporate the boundary conditions by expanding the fields as 

(/)a(x, t) =22 0n,«,aXn(x) (2.56) 

n,n ^ '^^ 

where aj„ == 7r(2n)/?i/3 for bosons and a;„ = 7r(2n + l)/fij3 for fermions. 
These are known as the even and odd Matsubara frequencies, respectively. 
Using this expansion we have 



(27ri)(i±i)/2 (n/3)J 



n,n.Q 



if we also take account of the jacobian involved in the change of integration 
variables. Note that the difference between the jacobians in the bosonic 
and fermionic case, is a consequence of the fact that for Grassmann vari- 
ables we have that J d(pf2(l) = ,(2 = /'^(/20)./2(/20) instead of the result 
J d(f2<j>)(l/ f2){f2'P) that we expect on the basis of ordinary complex inte- 
gration. Note also that these are all again gaussian integrals, so we find in 
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first instance 

^0 = n iPi-^^^n + en,a " A^))^' (2.58) 



exp < T X! ln(/3(-i?ia;„ + en,Q - A*)) 

l^ n,n,Q ^ 



To evaluate the sum over Matsubara frequencies we need to add a con- 
vergence factor e*'^"'' and finally take the limit ry | 0. The precise reason 
for this particular procedure cannot be fully understood at this point but 



is explained in section 2.5. However, doing so we indeed find that 



limyin(/3(-i?iw„ + e-^))e*"""=ln(lTe~'^('"''^) . (2.59) 

n 

To see that this is correct, we differentiate the latter equation with respect 
to /3^. This gives us 

1 ■r-^ e*'^"'' 1 

S m V ^'^n - (e - M)/?i " ^^P^^^^^T^ ' ^^■^'^^ 

which can be proved by contour integration in the following way. The 
function ?i/3/(e^^^=F 1) has poles at the even and odd Matsubara frequencies 
with residu ±1. Hence, by Cauchy's theorem the left-hand side is equal to 

1 /■ , e''^ ±1 

lim / dz 



>)j.o 27ri J (J z - (e - ^)/?i e^P^ ^ 1 

with C a contour that fully encloses the imaginary axis in the direction 
shown in figure 3. Adding to the contour C the contour C , which gives no 
contribution to the integral, we can again apply Cauchy's theorem to obtain 
the desired result. Why does the integration over the contour C vanish? 
The reason for this is that the integrant behaves as ±e~'-'''^^'''^°^^''/|2:| for 
Re(z) -^ oo and as — e^'^^^^-'/kl f^^' Re(z) -^ — oo. The integrant thus always 
vanishes much faster than l/\z\ on the contour C" for any < ry < fi./3. 

2.4.3 Green's function method 

The third and last method is simplest and, therefore, most common in 
practice. We first rewrite the partition function as 



Zo = j d[(i>*]d[4>] expi y / dr /dx (2.61) 

^ T. r '^^' J '^''' C(x,T)G-l,,(x,r;x',T')</)a'(x',r')l 

a' ^ ) 
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Fig. 3. Illustration of the contour integration that is required to perform the 
summation over the Matsubara frequencies. The black dots indicate the position 
of the poles in the integrant. 



and see G ^ as a 'matrix' both in spin space as in coordinate space. We 
then know that this gaussian integral is just 



Zo = [det(-G-i)]Ti = exp{TTr[ln(-G-i)]} 
Clearly, we have from the action Sq [(j)* , cj}] that 



(2.62) 



G-_i„,(x,r;x',r') 

1(8 n^v^ 



(2.63) 



h \ dr 2m 
or equivalently that 



+ y^"(x) + e„ - M ^ <5(x - x')<5(t - t')S^ 



dr 2m 



+ y^-(x)+e„-M G„,„,(x,r;x',r') 



(2.64) 



= -M(X - x')<5(t - T')<5a, 



which means that Ga.Q'(x, t;x',t') is a Green's function. We come back 
to its physical meaning shortly. For now we only need to see that the last 
equation is solved by 



Ga,a'(x, T;X ,T ) 



(2.65) 






-ihiVn + en,a — M 



Xn(x)Xn(x')- 



^ — iLOnir — r') 
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thus again 



Zo = cxp S T ^ ln(/3(-i?iw„ + £„,„ - n))> 

y iLi,n.a ) 



(2.66) 



exp<^T2^1n(lTe-'^^'^"-°-''0 



In principle, we have slightly cheated in the last step of this derivation, 
because equation (2.65|) shows that 



-tfiLUn + en,a " A* 

The problem with this last way of calculating the functional integral is 
that it does not correctly account for the jacobian involved in the change of 
variables from imaginary time r to the Matsubara frequencies cun- However, 
this problem can be avoided by calculating never a single determinant, but 
always the ratio of two determinants. In this manner the effect of the 
jacobian is exactly cancelled. We will see some examples of how this works 
in sections 2.6 and 2.8.3| . 



We now return to the physical meaning of the above Green's function. 
To see its meaning, we first consider the time-ordered expectation value 

(r[V'a(x,r)^l,(x',T')]> (2.68) 

= 0{t - r')(V^a(x, t)^PI (x', t')) ± e{T' - r){^Pl (x', r')^a(x, r)) . 

Here the expectation value is taken in the grand-canonical ensemble and 
■i/'Q,(x, r) is the imaginary time Heisenberg operator, which is defined by 
^{H-^iN)T/h^^^^^^-{H-fj.N)T/h g^j-^^ therefore obeys the Heisenberg equation 
of motion hdripai^, ''') = [H ~ fJ-N, ipai'x., t)]_ . For the noninteracting case 
it reads 



h—M^,r)^ i-^—-V''^{^)-e^ + fi\^^ix,T) . (2.69) 

As a result 

n^ {T[M^, t)^1 (x', r')]) = M(t - r')([^«(x, r), V^l, (x', r')]^) (2.70) 



(^-^"^(x)-6.+/.) (T[V:.(x,r)^t,(x',r')]) • 



Substituting the equal-time relations that we derived in section 2.1 

['(/'q(x, r), ^/j^, (x',r)]ip = (5(x — x')(5q,q', thus suggests that 



I.e., 



G„,a'(x,r;x',T') = -(T[^„(x,T)^^,(x',r')]) • (2.71) 
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We can actually also prove this important relation, that bridges the gap 
between the functional formulation of quantum field theory used here and 
the more familiar operator formalism. First of all, it is clear from the slicing 



procedure used in our derivation of the functional integral in section p.3| , 
that it automatically leads to time-ordered expectation values. We should, 
therefore, only be able to prove that 

-G<,,„-(x,T;x',r') - (0„(x, r)C,(x', r')) (2.72) 

^ 1-J d[<t>*]d[<t>\ 0„(x,T)C'(x',r')e-^«[**'*l/'' . 

This is most easily achieved in the following way. We introduce a partition 
function in the presence of external currents J a (x, r) and J* (x, t) , where 
in the fcrmionic case these currents are also Grassmann variables. The 
partition function is taken to be 

Zo[J.J*] ^ Jd[r]d[4,] exp J -is'o[0*,0] (2.73) 

+ 51/ ^■^ / cfx ((?!);(x,r)J„(x,T) + J*(x,t)0„(x,t)) I , 
since then we must simply prove that 



trQ,Q'(x, t;x ,T ) 



±1 S^Zo 



(2.74) 



Zo SJ^{^,T)SJa.'i^',T') 

Using a short-hand notation we have 
Zo[J,J*]= I dmd[(t>] e^^{{cj,\G-^\(l>) + {<j,\J) + {J\4>)} . (2.75) 



The terms in the exponent can be rewritten as (0 -I- JG\G^^\4) + GJ) — 
(J|G| J), which is usually called completing the square. Performing a shift 
in the integration variables, we then easily see that 

Zo[J, J*] = Zo[0,0]e"("'l'^l'^) = Zo[0,0]expi-^ / dr /dx (2.76) 

X ^ / rfr' / dx' J*(x, t)Gq,q'(x, T;x',r')Ja'(x',r') > , 

which after differentiation indeed leads to the desired result. We can now 
in fact calculate the expectation value of the time-ordered product of any 
number of operators. With the above expression for Zq[J, J*] one can easily 
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prove that this results in the sum of all possible products of time-ordered 
expectation values of two operators. For instance 



6:(X, T)CJ>1, (x', T')0a" (x", r")0a'" (x'", t'")) (2.77) 

= (C(x, r)0„., (x'", t"'))(C' (x', T')0a" (x", r")) 
±(C(x, t)0„. (x", t")) (C' (x', r')0„"' (x'", t'")) . 



This is the famous Wick's theorem, which plays a crucial role in the next 
section where we start to discuss the profound effects that interactions can 
have on the results obtained thusfar. 



2.5 Interactions and Feynmann diagrams 

The Green's function, or one-particle propagator, Gq.o'Cx, t; x',t') is one 
of the most important quantities that we want to determine theoretically, 
because it gives us the possibility to calculate in principle the expectation 
value of any one-particle observable. It also gives us the elementary excita- 
tions of our system of interest. For example, the average density of the spin 
state I a) is in an ideal gas given by 

(V't(x,T)^„(x,r)) =TG„,c«(x,T;x,r+) (2.78) 






2 



n,n 



where we introduced the notation r+ for the limit rj J, of r + 77. Note 
that this procedure is necessary due to the timc-ordcning involved in the 
definition of the Green's function and the fact that the field operators do not 
commute at equal times. It also gives a natural explanation for our previous 
use of the convergence factor e*'^"''. Moreover, from the argument of the 
Bose or Fermi distribution function, we see that the elementary excitations 
have an energy of £„,« — M- If we replace in the Fourier transform of the 
Green's function, i.e., in 

G'a,a'(x,x';iw„) = -Ma,a' V"^ ■ Xn(x)Xn(x') , (2.79) 

^ -iTlLJn + en.a - M 

iujn by Lo we also see that G'q_q,'(x,x';u;) has a pole at huj — e„Q, — /i. 
This is a general result: Poles in GQ,_a'(x, x';a;) correspond to the energies 
of the elementary excitations. These energies in general can also have a 
negative imaginary components, which then correspond to the lifetime of 
the excitation. The question that arises, therefore, is how to determine also 
the Green's function for an interacting system. 
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Fig. 4. Diagrammatic representation of the partition function up to first order in 
the interaction. 



This can be done in perturbation theory, as follows. We want to calculate 
the expectation value 



6,(x,T)C'(x',r')) 



(2.80) 



with the action 5 [(/)*, (/)] = Sq[(J)*,(J)] + Sint[4i* , (/)]■ We now expand both the 
numerator and the denominator in powers of S'int- Up to first order we find 
for the partition function 



Z 






(2.81) 



= ■Z'o ( 1 - T- (S'int [0*, '/'Do 



Using Wick's theorem, we thus have 
-\{SintWA])Q = \Y.j' dTfd^fd^'Go;c.A^,T;x,T+) (2.82) 

a. a' 

., ~nx-x'), 



-Go;Q',a'(x',T;x',r+) 



/•hp 



±-^ / dr dx dx' Go;a,a(x', r;x,T+) 

~V(x — x') 

X Go;a,a(x,r;x',r+) . 

n 

To understand the general structure of the perturbation expansion, it is 
very convenient to represent equations (2.81) and ( 2.82| ) in terms of Feyn- 
man diagrams. The final result is shown in figure 4, where a wiggly line 
corresponds to the factor —V{x — x')/Ti and a thin arrowed line point- 
ing from (x',r') to (x, t) represents the noninteracting Green's function 
G'o;a,Q(x, r; x', r'). In figure 4 we have, for clarity sake, also explicitly indi- 
cated the various coordinates and spin degrees of freedom that we have to 
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7 « 




X,X 





Fig. 5. Diagrammatic representation of tlie numerator of equation ( p.SOD up to 
first order in the interaction. 



either integrate or sum over, respectively. It is, however, much more com- 
mon in practice to suppress in Feynman diagrams those degrees of freedom 
that have to be integrated or summed over, and denote only the degrees of 
freedom on which the quantity of in terest depends explicitly. For example, 
for the numerator of equation (2.80) we obtain 



- / dmd[cb] 0„(x,T)C'(x',T')e-^°[*^^l/'' 1- 



Su 



(2.83) 



Zo { -(0a(x,T)C'(x',r'))o + ^(0a(x,r)^;,(x',T')5int[0*,</'])o 



whose diagrammatic equivalent is then given in figure 5. Up to first order the 
interacting Green's function, which is represented by a thick arrowed line, 
obeys the equation shown in figure 6. Note that the so-called disconnected 
diagrams have exactly cancelled. This is a general feature that happens 
for any expectation value, i.e., the disconnected diagrams that occur in the 
numerator are always exactly cancelled by the denominator |18|,|l9|]. 

One can also show that the general structure of the perturbation ex- 
pansion of the Green's function is as in figure 7, in which only one-particle 
irreducible diagrams contribute to the selfenergy h'S. These one- particle 



X, t X, T 

ail = a 

x'.t:' x',x' 





X ,T 



X ,T 



Fig. 6. Interacting Green's function up to first order in the interaction. 
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Fig. 7. Exact Dyson equation for the interacting Green's function. 

irreducible diagrams can be distinguished from the fact that they do not 
become disconnected if we cut a single thin arrowed line. This is Dyson's 
equation |2^. It is particularly insightful for a homogeneous gas or for a 
trapped gas in the weak-coupling limit, since then the one-particle states 
Xn(x) are not affected by the interactions and the exact one-particle prop- 
agator can be written as 

g-iw„(r-r') 

Ga,a'{:X^,T;x',T') = 2^Ga,a'{n,iuJn)Xn{x)Xn{'X.') ^ ■ (2.84) 

n,n 

Note that for realistic trapped atomic gases this weak-coupling limit is es- 
sentially always realized in the normal phase of the gas. It is, therefore, also 
of some experimental interest to consider that limit in more detail first. 
One then finds that the Dyson equation in first instance becomes 

Ga,a' (n, iuJn) = Go;Q,a' (n, iujn) (2.85) 

+ ^ Go-a,a"{n,iuJn)T.a"^a"'{n,iU!n)Ga"',a'{'n.,iuJn) ■ 
a" ,ct"' 

Because G'o;a,Q'(n, iwn) = G'o.jj(n, «ti;„)(5Q,_Q,/ and because we have assumed 
the interactions to be spin independent, we can easily convince ourselves 
that in perturbation theory it is always true that the selfenergy obeys 
SQ",Q-"'(n, iwn) = So,//(n, jtt'n,)(5a"^Q,'" and as a result also that up to all 
orders in the interaction G^.a' (n,za;„) — Ga{'n.,iuJn)Sa,a'- As we will see 
in section ^^, however, nonperturbative effect can change this result due 
to the phenomena of spontaneous broken symmetry. Ignoring the latter for 
the moment, we find for each spin state the uncoupled equation 

Ga{n,iujn) = Go;a(n, iw„) -h Go;a(n,iw„)S„(n,jw„)Ga(n,iw„) , (2.86) 

or 

^ /. ^ - ^ / . ^ - S,(n, iujn) (2.87) 

Ga[n,iuJn) Go;Q(n, zw„) 

n 
Hence 

Gc.(n, *^„) = — ~^ . , (2.88) 
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Fig. 8. Selfenergy up to first order in the interaction. 

which shows that the poles in the Green's function are indeed shifted by 
the interactions. 

In our lowest order calculation we have found for the selfenergy the dia- 
grammatic result shown in figure 8. This is the most simple approximation 
we can think of. In the weak-coupling limit it is mathematically equal to 

hT.a{n,iuJn) = y^ V„,n';n,n'^^7; — 1 (2-89) 



/ ^ v*'n,n';n,n' 
w'ol' 



\ / J •'n',n;n,n' 

n' 

^ *^n' ,n;n,n'^Q,a' j 



1 



1 



The first term is known as the direct or Hartree, and the second term as 



the exchange or Fock contribution to the selfenergy. From equation (2.89) 
we conclude that the matrix elements of the interaction enter only in the 
combination Vn,n';n,n' ± Ki',n;n,ii'<5Q;,Q'. Clearly, this is a refection of the 
Pauli principle, which forces the effective interaction between two fermionic 
atoms in the same state |n, a) to vanish. 

To make the theory selfconsistent we should use in the expression for 
the selfenergy not the noninteracting propegators, but precisely those that 
follow from the Dyson equation. This is the Hartree-Fock approximation, 
which is diagrammatically given in figure 9. It leads in the weak-coupling 
limit to the new dispersion relation hu!n,a = Cn.a + ^^^^(n, 0) — /i for the 
'dressed' particles or quasiparticles of th e gas, where the selfenergy is es- 
sentially also found from equation ( ^.89| ). We only have to replace in the 
right-hand side en', a' ~M by £n',a' +^5]^,^(n', 0) — /i. In this manner we have 




C V-^.^---^-^ + '^...---o 



Fig. 9. Hartree-Fock approximation for the selfenergy. 
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thus obtained an approximation to the interacting Green's function that is 
nonperturbative in the interaction and effectively sums an infinite number 
of Feynman diagrams. Of course, the Hartree-Fock approximation can also 
be used in the strong-coupling limit. In that case it diagrammatically still 
corresponds to the solution of the Dyson equation of figure 7, with a selfen- 
ergy as shown in figure 9. However, we are then no longer allowed to assume 
that the exact Green's function is diagonal in the eigenstates Xn(x) of the 
external trapping potential. To see more explicitly what this means physi- 
cally, it is convenient to rederive the Hartree-Fock theory without making 
use of perturbation theory. 

2.6 Hartree-Fock theory for an atomic Fermi gas 

Because the Hartree-Fock theory is very useful in a large number of circum- 
stances, and because we want to illustrate a very useful technique which 
is nowadays often used in the literature, we are now going to reproduce 
the above results without making use of a diagramatic expansion. We con- 
sider for simplicity a fermionic mixture with an equal number of atoms in 
two hyperfine states and start by splitting our spin-independent interaction 
T^(x — x') into two spin-dependent parts such that one contributes only to 
the Hartree diagram and the other only to the Fock diagram. Denoting a 
spin-dependent interaction by Va\i3i-^a,p{'^ — x') we thus want that 

y(x - x')<5.,a'<5^,/3' - T/„V;a,Mx - x') + Fi^,0.„,^(x - x') (2.90) 

with Y.(3 y"a-aA'^ - x') = E/3 ^a^/3;a,/3(x - x') = 0. Usiug Operators in 
spin space, a possible solution to these equations is 

f« = ^(2-A2)T^(x-x') (2.91) 

and 

l/^ = i(2A2-l)V^(x-x') , (2.92) 

where P12 = (1 + <ti • <T2)/2 is the exchange operator of the spins of atoms 
1 and 2, and cra,a' are the usual Pauli matrices. Clearly, we now have 

Sir,tWA] = \ Y. I drfd^fd^' (2.93) 



Q,a';/3,/3' 







{C' (x, T)(f>c (x, t)FJ_^,.„^^(x - x)(l)p, (x', t)4>p (x', r) 



which we write as 



5int[</'*,0] ^ ^ir^\v"\rci,) - ^ir^\\v^\\r<P) • (2.94) 
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We next apply a so-called Hubbard-Stratonovich transformation |^,|2^ to 
both the Hartree and the Fock parts of the interaction. First the Hartree 
part. 

We note that e"'^'"*^'^ '"^1 can be written as a functional integral over the 
four real fields contained in KQ,a'(x, r) = Ko{x,T)6a,a' + '^(x, r) • (Ta,a' as 



exp 



^(0*0|F^|(/.*0)1 (2.95) 



2h' 



/dMexp{^(.|y--|.)-^ 



i^\r'i^)-^{r(b\K) 



If we put S^^[(t)* , (j)] ~ for a moment, we can substitute this equality in 
the partition function to obtain Z = J d[K]d[(j)*]d[(l3]e-^^''^'^'^'^'^/'' with the 
action 



a, a' 



S{k,^\^\=--[k\V'' \n) + Y,J dry dxC(x,T)|U— (2.96) 



2m 



+ y''(x) + Eq. - ^ 6a,a' + Ha,a' (x, t) > (j)a' (x, t) 



We see that in this manner the action for the fermions has become quadratic 
with a selfenergy ?iI]Q^Q/(x, r; x', r') = /{^^^'(x, t)J(x — x')5(r — r'). Hence, 
we can now integrate out the fermion fields to obtain Z ~ j' d[K\e~^ ['^'/^ 
and 

^"''M = -li^V^'^t^) - aTr[ln(-G-i)] , (2.97) 

with G^^ — Gq — S also a functional of Kq.q'(x, r). Up to now we have 
not made any approximations and have only performed an exact rewriting 
of the partition function. However, the action 5°*^[k] contains all powers of 
the fields Kq,^q'(x, t) and is thus quite complicated. It may therefore appear 
that we have not made any progress. 

This is not the case, because the partition function will be largest for 
configurations that minimize the action 5°*^[k]. To make use of this obser- 
vation, we expand the action around its minimum, i.e., we put Ka,a'(x, r) = 
(kq,^c('(x)) -|-5kq,q'(x, t) and require that (55''^[k]/5kq.^q'(x, t)!^^^^) = 0. If 
we neglect the fluctuations, we find that Z ~ e^'^ [<«>]/''. This will turn 
out to be the Hartree approximation, as we will show now. First of all we 
have 

G"i ^ G^^ - {k) /fi- 5Klh= G"'^ - Sk/U ^ G"'^ il-G" Sk/K) . (2.98) 



Substituting this in equation ( ^.97 ) gives us for the terms linear in the 
fluctuations 

- nTv[-G"5K/n] - (5k\V"'^\{k)) . 
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If (Ka,Q'(x)) is indeed a minimum of the action S' [k], these hnear terms 
have to vanish, which imphes that 

(«:„,,„(x)> =^ /'dxt/i?,/3';„,M^-^')G^/3'(x',T;x',T+) . (2.99) 

As promised, this is precisely the most general mathematical expression 
for the Hartree contribution to the selfenergy in figure 9. Taking again 
('ta',a(x)) = k(x)(5q,'^q,, we also have that the Hartree approximation to the 
one-particle propagator obeys G^^,(x', t; x', r+) = n(x')(5^^^//2, with n(x') 
the average total atomic density, and we obtain the usual Hartree selfenergy 

.c(x) =^ /dx' Fi^^^^^^(x-x')^ = /dx' y(x-x')n(x') . (2.100) 

Now we again want to include the Fock part of the interaction and 
treat this also by a Hubbard-Stratonovich transformation. This requires 
introducing four real fields that depend on two spatial coordinates. They 
are denoted by Aa.Q'(x,x',T). Schematically we then use 

exp|^(r0||F^||r(/.)| (2.101) 

d[X] exp |- JL(A||T/^-^||A) + ^(A||r0) + ^irm) 

which leads to the total selfenergy 

nS„,„,(x,T;x',T') = K„,(x,T)(5(x-x')-A„^„,(x,x',T)]r5(r-T') (2.102) 
and after integration over the fermion fields to the effective action 

S''"[>.,\]=-\{K\V"-'\^) + ^iX\\V^-'\\X)-hTr[H-G-')] . (2.103) 

Requiring now that also 6S'^^[K,X\/SXa,a'{x,^' ,t)\\^^\) = 0, then leads to 
the expected Fock selfenergy 

(A„.,„(x',x)) - 5] /dx' y^^,„,.„^^(x-x')G^,S(^''^;^'^^) • (2-104) 

In the symmetric case it simply becomes (Aa'^Q(x',x)) = X{x' ,x.)Sa'.a with 

A(x',x) = T/(x-x')n(x',x) (2.105) 

and n(x',x) ^ J:^GI${^' ,t;^,t+)/2 = E^(V'^(x, r)V'^(x', t))/2 exactly 
the Hartree- Fock approximation to the off-diagonal part of the one-particle 
density matrix. 
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Note that to actually perform the Hartree-Fock calculation in the latter 
case, we need to be able to determine the Green's function G^^, (x, r; x', r'). 
The easiest way to do so, is by realizing that it is the Green's function of 
the operator in the fermionic piece of the action S[k, A, (j)* , (j)] obtained after 
the Hubbard- Stratonovich transformations. If we diagonalize this operator 
by solving the eigenvalue problem 



+ y-(x)+.j(x)-6'„ x'n(x) (2.106) 







2m 

-y'dx'A(x,x')xUx') 

the desired one-particle propagator aquires the ideal gas form 

Gf_f,(x,T;x',T') (2.107) 



^a,a' / ^ 



-xn(x)x::(x'^ 



-iflLOn + £„ a - M " " ^/3 



with new one-particle energies e^ „ = e^ -I- Cq and eigenstates Xn(x) that 
incorporate the average effect of the interactions of an atom with all the 
other atoms in the gas. Clearly, in the cases that the eigenstates are not 
affected by these so-called mean-field effects, we recover the weak-coupling 
results of the previous section. 

Although we have thus precisely reproduced our diagrammatic result, 
there are two important advantages in using the Hubbard-Stratonovich 
transformation. First, it is in principle exact, and allows us to also cal- 
culate corrections to the Hartree-Fock approximation. For example, if 
we expand S" [k. A] up to quadratic order in 5k and (5A, and neglect all 
higher orders, we find the so-called Generalized Random Phase Approx- 
imation. The latter approach actually gives us also the opportunity to 
study the density fluctuations and therefore the collective excitations of the 



gas. We return to this important topic in section 3^. Second, it allows 
for a beautiful way to describe phase transitions. As mentioned previ- 
ously, in perturbation theory we always find for a homogeneous gas that 
G'Q,a'(x, t; X, r+) = n5a,a' /2 due to the same feature of Go;a,Q'(x, r;x, t"''). 
From a fundamental point of view this is a result of the translational in- 
variance of the gas and of the rotational symmetry in spin space. However, 
we can imagine that in principle we can also find sclfconsistent solutions 
that do not have this property. We then have a spontaneous breaking 
of symmetry and therefore a phase transition in our system. For exam- 
ple, if below a certain temperature Ga.a'{'^,T]ii.^T^) — n('x.)6a,a' /'^, we 
are dealing with a transition to a charge density wave or (Wigner) crys- 
tal. If on the other hand Gq,,q'(x, t;x,t"'") = ?t.5q_q,'/2 -|- m • cTq, „/ the gas 
is in a ferromagnetic phase. For a spin-density wave we even have that 
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G'q.q'(x, t; X, r+) = n(5Q,,a//2 + m(x) ■ CTa,a' ■ In all these cases the Hubbard- 
Stratonovich approach used above leads in a natural way to the appropriate 
Landau theory of the phase transition. Since the Landau theory is also 
very useful for the understanding of superfluidity in atomic gases, which is 
clearly the primary goal of the present course, we use the next section to 
give a short introduction into this subject. 

2.7 Landau theory of phase transitions 

In this summary of the Landau theory of phase transitions, we restrict our- 
selves in first instance to the homogeneous case because this makes the 
discussion more transparent. However, at the end of the section we also 
briefly mention how the inhomogeneity enters the theory. We have seen, for 
example in the Hartree-Fock theory discussed above, that the Green's func- 
tion Ga^a'(x, t; X, r+) ~ nda,a'/'2 + (s) • cTa^a' Can have a nonzero value of 
the average spin density (s) — {ipl^{x,T)cra.a'ipa'{^,'''))/2, which is usually 
called the magnetization m in this context. This signals a phase transition 
to a ferromagnetic phase and the magnetization is called the order param- 
eter of this transition. In the previous section we have also seen how we 
can, in principle exactly, obtain an expression for the partition function as 
a functional integral over the field s(x, t), i.e., 

d[s] e-^°"H/?^ . (2.108) 

We just have to integrate out the fields kq and Xa.a'- For space and time 
independent values of the magnetization s the action 5''^*[s] is h(3VfL{s), 
where /l(s) is the Landau free-energy density and V is the total volume 
of system. Because of the symmetry of S[4>*,4'] under rotations of the 
spin, this free-energy density must also be rotationally invariant and can 
therefore only depend on the magnitude of s, which we denote by s. If a 
phase transition occurs the behaviour of /l(s) can essentially fall only into 
two categories. 

At temperatures very high compared to the critical temperature Tc, the 
system is fully disorded and the free energy /l(s) must have a single mini- 
mum at s =: to make sure that the order parameter (s) is zero. Bringing 
the temperature closer to Tc, however, the free energy can develop a second 
local minimum. As long as the free energy in this second local minimum 
is higher than the minimum at s — 0, the equilibrium value of (s) will still 
be zero and no phase transition has occured. Lowering the temperature 
further, the value of the free energy in the second minimum decrease until, 
precisely at the critical temperature T^ it is equal to the free energy at 
s = 0. For temperatures below this critical one the second minimum has 
actually become the global minimum of the free energy, which implies that 
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Fig. 10. Qualitative behaviour of a) the Landau free energy and b) the order 
parameter for a discontinuous phase transition. 



(s) 7^ and we are in the ordered phase. In this scenario the order pa- 
rameter has always a discontinuity at the critical temperature. As a result, 
this corresponds to a discontinuous, or first-order, phase transition. The 
behaviour of the free energy and the order parameter is illustarted in figure 
10. This should be compared with the behaviour of the free energy and the 
order parameter for a continuous, i.e., second or higher order, phase transi- 
tion, which is quite different and depicted in figure 11. Now the Landau free 
energy fL{s) has always a single minimum, which for temperatures above 
Tc is at s = but at temperatures below the critical temperature shifts 
to a nonzero value of s. In particular, the order parameter thus shows no 
discontinuity at T^. 

In the case of a second-order phase transition, wc conclude that near the 
critical temperature (s) is very small. As a result we can then expand the 
Landau free-energy density in powers of s. Because the free energy must 
also be symmetric under s ^ — s, we have 



/i(s) = «(T)|s|2 + ^|s| 



(2.109) 




with /3 > and a{T) = ao{T/Tc - 1). Thus if T > Tc we have a(T) > 
and (s) = 0. But for T < T^ (x{T) becomes negative and we have 



(2.110) 



Note that the free-energy density in this minimum is — q;q(1 — T/rc)^/2/3 < 
and has a discontinuity in its second derivative with respect to the tempera- 
ture. Historically, this is the reason why the corresponding phase transition 
was named to be of second order. 

In fact, the Landau theory of second order phase transitions is slightly 
more involved, because it also considers slow spatial fluctuations in s. Since 
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Fig. 11. Qualitative behaviour of a) the Landau free energy and b) the order 
parameter for a continuous phase transition. 



a uniform rotation of s costs no energy, we must have that the free energy 
is 

FM= j d^ (^7(Vs(x)f + a(T)|s(x)|2 + ^|s(x)|4) (2.111) 

and Z ~ /cf[s] e^^^'^^^\ Landau theory now escntially corresponds to min- 
imizing Fl[s\ and taking only quadratic fluctuations into account. Taking 
also higher order fluctuations into account turns out to be very difflcult. 
It requires Renormalization Group methods, which we are not going to 
treat here and that can be found in the literature |jlq] . The effects of these 
fluctuation corrections to Landau theory are generally known as critical 
phenomena. Interestingly, they do not occur in first-order phase transi- 
tions. More importantly for our purposes is, however, that for atomic gases 
the fiuctuations are only important in a small temperature interval around 
the critical temperature. For many applications it is, therefore, possible to 
neglect them. What usually cannot be negected is the effect of the inhomo- 
geneity of the gas. In the context of Landau theory that implies that the 
coefficients a{T), (3 and 7 in equation ( 2.111 ) become also dependend on 
the spatial position in the trapping potential. We will see several examples 
of this feature in the following. 



2.8 Superfluidity and superconductivity 

We finally want to consider two important second-order phase transitions, 
that are purely due to quantum effects and can be conveniently treated with 
the methods that we have developed sofar. Moreover, they occur very often 
in nature, for example in metals, in liquid helium and recently of course also 
in atomic gases of rubidium [no| , lithium ll|] , sodium |l3] and hydrogen [E3l . 
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2.8.1 Superfluidity 

Let us first consider the last case, which is associated with a gas of spin-less 
bosons. For the low temperatures of interest the action is 

5[0*>] (2.112) 

= y^ dry dx0*(x,T)|?i— -^— + y-(x)-A^J0(x,T) 
1 /■'"^ /■ 

+ - / dT dx Vo(f>*{x,T)(j>*{x,T)(P{x,T)4){x,T) , 

where V(x— x') is taken to be equal to Vb(5(x— x'). The justification for this 
simplificaties is, roughly speaking, that the thermal de Broglie wavelength 
Ath = (27r7i /mfcsT)^/^ of the atoms is much larger than the range of the 
interatomic interaction. In this system the phase transition of interest is 
Bose- Einstein condensation. The associated order parameter is (0(x, r)), 
since for time independent (/)(x, t) the above action has presicely the form 
of a Landau theory with a 'free-energy' of 

Fl[(I)*,(J)\ (2.113) 



Jdx (^|l|V0(x)p + (y-(x) - m)|0(x)P + ^|0(x)| 



We conclude that in lowest order the critical temperature is determined by 
fJ,{Tc) = €o, because then the configurations (/)(x) ex Xo(x) precisely make 
a vanishing contribution to the quadratic part of the 'free energy'. This 
condition makes sense, because it is exactly the condition that we have 
a Bose-Einstein condensation in the ideal case: In the ideal Bose gas the 
number of particles in the one-particle ground state is A^o = l/(e^^'^°~'^'' — 1), 
which indeed diverges for /x = eq. 

To determine the corrections to this result we now explicitly substitute 
0(x, r) = 0o(x) -I- 0'(x, t) into our functional integral. It is for lateron 
important to realize that to consistently define the fluctuations (/)'(x, r) in 
this manner, we also have to require that 

/dx 0S(x)0'(x,r) + f dx (^o(x)0'*(x,t) = . (2.114) 

The physical reason behind this condition is that (p'ix, t) should contain all 
the configurations that are orthogonal to 0o(x). In principle, therefore, it 
should also contain configurations that in effect only multiply (/>o(x) by a 
global phase. Such fluctuations lead to the phenomena of phase 'diffusion' 



and are discussed in section |3.2| . In full detail we find after the above 
substitution that 

5[</)'*,0']=n/3Fi[0S,0o]-f 5*0 [./)'*, 0']+5i„t[0'*,</''] , (2.115) 
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where the Unear and quadratic terms are given by 
SoW*A'] (2.116) 

= J drjd^ cl>'*{^, t) \-^^ + V^°"(x) ~fi + yo|(^o(x)P Uo(x) 

2m 



+ £ drjd^ 0'(x,t) |-^ + ^^'^(x) -^i + Vo\M^)\^ [ '^S(x) 

+ dT dX(f>'*{x,T) 

X l^i- - ^ + ^'"W - M + 2Vo\M^)\' \ 0'(x,r) 
1^ OT 2m ' 

1 r^^ f 

+ -J drj dx Vo(0o(x))20'*(x,t)0'*(x,t) 



1 



R/3 



+ ^J drJ dx yo(0S(x))'0'(x,T)0'(x,T) , 
and the cubic and quartic terms by 
SintW*,cb'] = / dr / dx Vb0o(x)(^'*(x,r)(^'*(x,r)(^'(x,T) (2.117) 

+ / dr / dx T/o0S(x),/)'*(x,r),/)'(x,T)0'(x,T) 

1 /■'''^ /■ 
+ -y drj dxyo0'*(x,r)0'*(x,r)0'(x,r)0'(x,T) , 

respectively. 

In the Bogoliubov approximation we neglect the last three interaction 
terms [Q. Furthermore, to make sure that (0(x, r)) = 0o(x) or (0'(x, r)) = 
0, we need to require that the terms linear in 0' and (p' drop out of the 
action Sq[cJ)'* , (/)']. Clearly, this implies that 

2m 



+ y'=''(x) + Vo\M^)r M^) = ^^M^) , (2.118) 



which is the same result as obtained from minimizing the Landau 'free- 
energy' Fl[(P*,(J)]. In the context of trapped atomic gases, this equation is 
known as the Gross-Pitaevskii equation [E5| . It determines the macroscopic 
wave function of the condensate. The reason for calling 0o(x) the macro- 
scopic wave function follows from the fact that the total density of the gas 
now obeys 

n(x) = (0(x,r)<^*(x,r+)) = \M^)\' + (<^'(x, r)0'*(x, r+)) . (2.119) 
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The total numb er of condensate atoms thus equals Nq = J dx \(J)q{'x.)\'^ . As 
equation (2.119) shows, it is in general always smaller that the total number 
of atoms in the gas due to the effect of the fluctuations. Note that in our 



present formulation the average ' 



(x, T+)) physically describes not 



only the depletion of the condensate due to the usual thermal fluctuations 
known from the ideal Bose gas, but also due to the interactions, i.e., purely 
due to quantum fluctuations. 

Assuming that we have solved for the Gross-Pitaevskii equation, the fluc- 
tuation corrections are in the Bogoliubov theory determined by a quadratic 
action of the form 



Sb[ 



(2.120) 



= ~J rfrydx [0'*(x,t),0'(x,t)].G-i 



(/.'(x,r) 
0'*(x,r) 



where the associated Green's function G has now a matrix structure because 
not only the normal average ((/)'(x, t)0'*(x',t')) but also the anomalous 
average (0'(x, r)(/)'(x', r')) is now unequal to zero. We thus have that 



G(x,t;x',t') = 



0'(x,t) 
0'*(x,r) 



[^'*(x',r'),0'(x',r')] 



(2.121) 



From equation (2.116) we in fact find that in the Bogoliubov approximation 
G-i(x, r; x', r') = G^\x, r; x', r') (2.122) 



Fo(0S(x))2 2Vo|0o(x)|2 
with the noninteracting Green's function Gg defined by 



(5(x - x')(5(t - r') 



Go-1(x,t;x',t') = 



Go\x,T;x',r') 

Go-i(x',T';x,r) 



(2.123) 



and, of course, 




Go-1(x,t;x',t') 




= nl^dr' 


2m 



(2.124) 



+ y°'=(x) - ^ U(x - x')<5(t - r') 



This is clearly only the lowest order result for the Green's function, because 
a perturbative treatment of Sint[(t>'* ,(/>'] leads, in the same way as in section 
2.5, to higher order corrections. In general, the Dyson equation is, however, 
always of the form 



Gil Gi2 
G21 G22 



-I -1 



Gn-^ 







\ Sii 


S12 





Go-^ 




S21 


S22 



(2.125) 
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Fig. 12. Exact Dyson equation for the interacting normal and anomalous Green's 
functions. 



The off-diagonal elements are again called anomalous, since they vanish in 
the normal phase of the gas. Diagrammatically the Dyson equations for 
Gil and G21 are shown in figure 12. 

The selfenergy in the Bogoliubov approximation is 



?iS(x, r;x', r 






6{x-x')6{t-' 



(2.126) 

(2.127) 



Because we have completely negelected the cubic and quartic terms in the 
action to obtain this selfenergy, the Bogoliubov approach is only valid if the 
fluctuations are sufficiently small. Physically, this impies that the deplec- 
tion of the condensate must be small. It can, therefore, not be applied to 
liquid helium, but is expected to be valid for a weakly-interacting atomic 
gas at such low temperature that it essentially only consists of a conden- 
sate. Under these conditions the Bogoliubov theory not only predicts the 
condensate density profile rio(x) — |(/)o(x)p but also the collective modes 
of the condensate. Both these predictions have been accurately verified in 
recent experiments p^-[2^] . Theoretically, the eigenfrequencies of the col- 
lective modes are again determined by the poles in the one-particle Green's 
function. To understand more clearly how these can be determined, we first 
consider a homogeneous Bose gas. 

In a box with volume V = L^ the one-particle states are most easily char- 
acterized by the wavevector k = (27r/L)n and equal to e^^'^/VV. The one- 
particle energies are thus eis, = h k^/2m. Moreover, the Gross-Pitaevskii 
equation reduces to ^ = Vb|0op in that case, because the Landau 'free- 
energy' is minimized for a macroscopic wave function that is independent 
of the position in the box. Making use of this fact, equation ( 2.122| ) can 
immediately be solved by a Fourier transformation. The result is 






(2.128) 



Vol 



ifwo. 



+ ek + V'o| 
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Clearly there are poles in G(k, uj) if the determinant of the right-hand side 
is zero or if 



hio = hcji, = Je^ + 2Vb|0oPek = \hi + 2Vonoek 



(2.129) 



This is the famous Bogoliubov dispersion of the collective excitations. Note 
that to finish the calculation we still have to obtain the condensate density 
no — \(j)o\'^- This is determined by the total density of the gas, which obeys 

n = |0oP - G'ii(x, r;x, T+) = no + n' , 



(2.130) 



with n' the density of the 'above' condensate particles. By inverting the 
right-hand side of equation (2.128), we find that the noncondensate density 
is equal to 



lim ■ 






j ihun + (^k + Vono 



(2.131) 



V ^ 

kT^O 



fk + V"ono 1 



fUjJ\, 



j/JRWk: _ \ 



Ck + Vpno - hcji^ 



For a given density and temperature, the last two equations thus fully de- 
termine the condensate density. Note that equation ( 2.131 ) explicitly shows 
that the condensate is indeed depleted both by thermal as well as quantum 
effects. 

The generalization to the inhomogeneous case is straightforward. First 
we again have to solve the Gross-Pitaevskii equation at a fixed chemical 
potential. Given the condensate wave function, we can then calculate the 
collective mode frequencies by finding the poles of G, or equivalently but 
more conveniently, the zero's of G~^. Clearly, the latter are located at 
huj = hujn, where htOn is found from the eigenvalue problem 



vom^)? 



Vo(0o(x))2 



Mn(x) 
'Un(x) 



(2.132) 



huJri 



-1 



M„(x) 
l'n(x) 



and we introduced the operator K — ~h V /2m + V^'^'^(x) — /i. This is 
de Bogoliubov-de Gennes equation that has recently been applied with 
great succes to the collective modes of a Bose condensed rubidium and 
sodium gas |29, B0[. Note that a special solution with hwo = is given by 
[uo{x),vo{x)] = [0o(x), — 0q(x)]. As a result we can, by making use of the 
fact that the left-hand side of the Bogoliubov-de Gennes equation involves 
a hermitian operator, easily proof that all the solutions obey 



dx (/>o(x)Wn(x) + dx 0o(x)Wn(x) = 



(2.133) 
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as required by the condition in equation (2.114). Moreover, we can similary 



show that the solutions with hui^ > can always be normalized as [ pl[ 

y"dx(|?/„(x)|2-|t,„(x)n = l. (2.134) 

Physically, the zero frequency solution n = describes the dynamics 
of the global phase of the condensate Q. Because of the so-called C/(l) 
symmetry of the action, i.e., its invariance under the simultaneous phase 
changes ^(x, t) -^ e*''^(x, t) and (/)*(x, r) -^ e~^^ (p* {x, t) , this solution is 
essentially of no importance to the thermodynamic properties of a macro- 
scopic gas sample and is therefore usually neglected. Nevertheless, it has 
from a fundamental point of view some interesting consequences, as we will 



see in section 3.2. Knowing all the eigenstates of G~^ we can then eas- 
ily perform the inversion and finally again determine the density profile 
of the noncondensed atoms. Keeping the physical significance of the zero 
frequency mode in mind, we ultimately find, 

n'(x) ^ Y. ((l"n(x)P + K(x)P) ^^J^_^ + b„(x)p] , (2.135) 

n#0 ^ ^ 



which may be compared with equation (2.131). 

At temperatures near absolute zero, we have as a good approximation 
that n'(x) = and the Bogoliubov theory applies. However at nonzero 
temperatures we thermally excite particles and n'(x) becomes nonzero. If 
we treat the effect of the noncondensate part of the gas in the Hartree-Fock 
approximation, wc find that 5'int [0' , 0'] on average adds 



S^J [(/.'*, (/.']- 2 y drj dx yon'(x),/)'* (x, r),/)o (x) 
+2 dr / dx Vbn'(x)0'(x,r)0;(x) 



(2.136) 



+2 / dT dx Von'{x)(f)'*{x,T)4)'{x,T) , 

to the action SqIcJ)' , (j)']. In figure 13 we indicate how this can be understood 
diagramatically. Performing the same analysis as before, we conclude that 
the Gross-Pitaevskii equation is modified to 

-^ + ^'="(x) + 2Von'(x) + K)|0o(x)P M^) = ^^M^) , (2.137) 

\ 2m / 

and the normal self energies are changed into 2Vb|(/'o(x)p + 2Von'(x) = 
2Vbn(x). The Bogoliubov-de Gennes equation for the elementary excita- 
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a) 



b) 
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Fig. 13. Hartree-Fock corrections to a) the hnear and b) the quadratic interaction 
terms of the Bogoliubov theory. This represents the Popov theory. 



tions is, therefore, now given by 
k + 2Vbn(x) 



1^(0O(X))2 1 




u„(x) 




fs: + 2yon(x) 




Wn(x) 




= huj^ 


■ 1 







L 






(x) ■ 



(2.138) 



These last two equations in combination with equation ( 2.135| ) are known 
as the Popov theory in the recent literature 133]. It is much studied at 
present in the context of Bosc-Einstein condensation in atomic gases, and 
has been applied with succes to the equilibrium density profile of the gas 



below the critical temperature |34|. It has also been used to determine the 
collective mode frequencies of the gas at nonzero temperatures, however, 
with much less success [§g,|3§] . The reason for the failure of the Popov theory 
in this case is that the Bogoliubov-de Gennes equation in equation ( ^.13§| ) 
describes physically only the motion of the condensate in the presence of a 
static noncondensed cloud and not the dynamics of the noncondensed cloud 
itself. How that can be achieved is discussed lateron in section 3.5, when we 
have obtained a better understanding of the nonequilibrium theory. Now 
we briefly want to make a connection to the interaction parameter Vq used 
in the Bogoliubov and Popov theories and the specific two-body scattering 
properties of the atomic gas of interest in a particular experiment. 



2.8.2 Some atomic physics 

In the previous section wc mentioned that the Bogoliubov theory, and to a 
certain extent also the Popov theory, agree very well with experiment. It 
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Fig. 14. T-matrix equation for the effective interatomic interaction. 



is clear, however, that to apply these theories to an actual experiment we 
need to know the interaction parameter Vq. In the case of atomic gases 
it is indeed possible to perform an ab initio calculation of this quantity, 
something which for instance cannot be done for liquid helium. The reason 
why Vq can be determined for atomic gases, is that under the experimental 
conditions of interest the densities are always so low that we only need to 
consider all two-body processes taking place in the gas and we can neglect 
three-body and higher-body processes. This implies physically that we only 
have to calculate and add the quantum mechanical amplitudes for two atoms 
to scatter of each other an arbitrary number of times. Diagrammatically 
the procedure can essentially be summarized by the T-matrix equation in 
figure 14, because by iteration of this equation we easily see that we are 
indeed summing all two-body interacion processes. 

Denoting the total momentum of the two incoming particles by KK. and 
the sum of the two Matsubara frequencies by f7„ , the T-matrix equation in 
figure 14 can, with our knowledge of how to perform the Matsubara sum in 
the right-hand side, easily be shown to be mathematically equivalent to 



Vn 



T{K,inr,) = Vo + ^Yl 



1 



V ^ ihQn - eK/2+k - eK/2-k + 2At 



(2.139) 



X 1 + 



,/3(<SK/2+k-A') _ I g/3(eK/2-k-M) _ I 



T{K,in, 



It is immediately solved by 
1 11^ 



1 



T{K,in, 



Vo V ^-^ ihn,, 

k 



X 1 + 



^ eK/2+k — eK/2-k + 2^ 

1 1 



(2.140) 



g/3(cK/2 + k-A') _ I g/3(<SK/2-k-M) _ I 
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Note that this is a somewhat formal result, because the momentum sum 
in the right-hand side has an ultra-violet divergence, which from a field- 
theoretical point of view requires a renormalization procedure. 

The physical reason for the divergence is that we have used the ap- 
proximate potential Vo5(x — x') instead of the actual interactomic potential 
y(x — x'). The argument for using the 5- function approximation was that 
the thermal de Broglie wavelength of the atoms is for the ultra-low tempera- 
tures of interest always much larger than the typical range of the interaction. 
We now see that this argument is not fully correct, because if we calculate 
corrections in perturbation theory, we have to deal with momentum sums 
which are not always restricted to momenta of order /i/Ath and arc therefore 
sensitive to the precise details of the interaction potential. 



To cure this disease we note that if we put in equation (2.139) the 
Bose occupation numbers equal to zero, we precisely get the T-matrix or 
Lippmann-Schwinger equation p7| ] for two atoms interacting with the po- 
tential Vo(5(x — x'). Furthermore, if the atoms interact with the potential 
y (x— x') the solution of the corresponding Lippmann-Schwinger equation is 
known from elementary scattering theory to be equal to inah jra for small 
incoming momenta and energies, where a is the s-wave scattering length. 
We thus conclude that we must interpret 

11.^ 1 



Vo ^ k *^^" ~ eK/2+k - eK/2-k + 2/x 



k 

lo 'in/'i'nar 
principle obeys 

1 m 

T(K, zfi„) 47ra?i^ V ^ iMln ~ eK/2+k " eK/2-k + 2^* 

1 1 



as being equal to m/A'KaTi and, therefore, that the desired T-matrix in 

IV ^- ('2 141) 



g/3('!K/2 + k-A') _ \ g/3(eK/2-k-M) _ \ 

To obtain an accurate theory for a trapped atomic Bose gas that in- 
cludes the effect of all two-body processes, we should now use instead of 
Vq the above T-matrix in the Bogoliubov or Popov theories. Of course, 
in the normal phase the same is true for the Hartree-Fock theory. If the 
temperature is not too close to the critical temperature for Bose-Einstein 
condensation, it turns out that T(K, if2„) ~ 47ra?i /m and we arrive at the 
conclusion that for an application to realistic atomic gases we must replace 



everywhere in section 2.8.1 the interaction parameter Vo by Anah /m. For 
temperatures close to the critical temperature, this is however no longer 
true and the Bose distribution functions have a nonnegligible effect of the 
T-matrix |38,39[. The theory that, by using 'dressed' one-particle propa- 
gators in figure 14, selfconsistently includes these effects of the medium on 
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the scattering properties of the atoms is known as the many-body T-matrix 
theory. It results in an effective interaction that, to avoid confusion, is de- 



noted by r(+)(K,r2) in section 3.4. As we will see then, it plays a crucial 
role if we want to arrive at a consistent discription of condensate formation 
in an atomic Bose gas. 

2.8.3 Superconductivity 

Finally, we want to briefly discuss the Bardeen-Cooper-Schrieffer or BCS 



theory of superconductivity 1 40 1 , which has received considerable attention 
recently in connection with ongoing experiments with the fermionic isotope 
of lithium. The reasons for this will become clear shortly. We are in that 
case dealing with effective spin 1/2 fermions and the action is taken to be 

SW,<f\ (2.142) 

= H y^ dry dx<^;(x,r)|?i— -^— -^|0„(x,t) 

+ dr dx Vo(l)* (x, r)<?!)J (x, t)0| (x, r)(/i| (x, r) . 

Note that the use of a single chemical potential implies that we only consider 
the optimal case of an equal density in each hypcrfine state. This situation 
is optimal in the sense that, for unequal densities in each hyperfine state, 
it is no longer possible to pair up all the atoms in the gas. As a result the 
critical temperature of the gas drops dramatically, and in general becomes 
experimentally inaccessible. The effect is to a large extent analogous to 
putting superfluid ^He in an homogeneous magnetic field. Furthermore, 
we in first instance consider the homogeneous case, because we want to 
illustrate in this section the local-density approximation to include the effect 
of the external trapping potential. Physically, this approximation treats the 
gas as consisting of a large number of independent gases that are in diffusive 
equilibrium with each other. Such an approach only works if the correlation 
length of the gas is much smaller than the typical length scale associated 
with changes of the external potential. Fortunately, this is almost always the 
case for realistic trapped atomic gases and for that reason the local-density 
approximation is often used in practice. 

The BCS theory is the theory of Bose-Einstein condensation of so-called 
Cooper pairs. This means that the order parameter is (0|(x, t)0|(x, r)), in 
analogy with the order parameter ((/>(x, r)) for the Bose case just discussed. 
Furthermore, it requires that the interaction parameter Vq is negative, since 
otherwise the formation of pairs would not be energetically favorable. From 
now on we assume, therefore, that this is the case. The condensate of 
Cooper pairs can also be nicely treated with a Hubbard-Stratonovich trans- 
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formation. We now introduce a complex field A(x, r) and use 



exp ■ 



up 



dT / dx Vor^ (x, T)ri (x, r)0x (x, T)0t (x, r) (2.143) 



hf3 



d[A*]d[A] exp<' i ^dr f dx f\A^^hl)l. 



+ A*(x,T)0|(x,r)</)|(x,T) +</)*(x,r)</)*(x,T)A(x,r) 

This leads to a partition function with the action 

r^^' r |A(x,r)|2 



S[A*,A,(j)*,(j)] =- dr dx 



Vn 



(2.144) 



-J dT j dx [0I(x,r),</)t(x,r)] -G-i 



'/'i(x,r) 
'/'*(x,r) 



where 



G-i(x,T;x',T') = Go-i(x,r;x',r') 



(2.145) 



A(x,t) 

A*(x,r) 



(5(x - x')5(t - r') , 



and the noninteracting Green's function Gq is defined by 



Go'(x,T;x',r') = 
and 

Go-i(x,r;x',r') = - 



Go^(x,t;x',t') 

-Goi(x',T';x,r) 



(2.146) 



Wl-"^-^"^'^--'^'^^-^'^- ^'-'''^ 



We thus see that the fermionic part has exactly the same matrix struc- 
ture as in the case of a condensed Bose gas, only the selfenergy is now 



(5(x - x')(5(t - t') . (2.148) 



If we again integrate out the fermion fields, we get the effective action 



11 


Sl2 


1 





A(x,T 


21 


S22 


^ h 


A*(x,t) 






S^«[A*,A] = - dT 



|A(X,T 



riTr[\n{-G-^)] , (2.149) 
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which we can expand in powers of A by using G"'^ = Gq^ — S = G,^^(l — 
GqS) and therefore 

CXD _. 

- hTi-[ln{-G-^)] = -hTT[\n{-Go^)] +TiJ2 — Tr[(GoS)™] • (2.150) 



m 

m—l 



ExpHcite calculation Q shows that for space and time independent A we 
obtain a 'free-energy' density of the form of the Landau-theory of second- 
order phase transitions, i.e., 



,2 . ,r^n^ 7C(3) 



f^(\A\)=a{T)\A\'+NiO)j^^\A\\ (2.151) 

with as expected 

a(r) = -— + -V— ^ r f 1 - ^, ^ , ^ (2.152) 



Vo V ^ 2(ek - m) V e/^l^k-p) + i 

^7V(0)ln(^|^)^7V(0)(^|^-l) ■ (2.153) 

Here A^(0) is the density of states of a single spin state at the Fermi energy 
e_F = ?i kp/2m and the critical temperature is given by 

T. = ^^|^cxp|--^|, (2.154) 

n Kb I 2kp\a\ J 



if we use the same renormalization procedure as in section 2.8.2| to elimi- 



nate the interaction parameter Vq in favor of the negative s-wave scattering 
length a. 

Below Tc we thus have a nonzero average (A(x, r)) = Aq. Using this in 
the Green's function for the fermions and neglecting fluctations, we find in 
momentum space that 



~hG-'{k,iLu,,) = 
and therefore poles in G(k, lj) if 



ihuJn + Ek - M ^0 

A2 -{ihujn + (k - ^J■) 



(2.155) 



hiUk - V(ek-M)' + |Aop . (2.156) 

This disperion relation has clearly a gap of magnitude |Ao| at the Fermi 
surface. As a consequence |Ao| is known as the BCS gap parameter. Note 
that the critical temperature below which the gap becomes nonzero depends 
exponentally on the parameter 7T/2kp\a\. For typical atomic gases this 
is a quantity which is much larger than one, and the BCS transition is 
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experimentally inaccesible. The only exception at the moment appears to 
be ^Li, with its anomalously large and negative triplet scattering length of 
—2160 ap p2[ . This explains the present interest in spin-polarized atomic 
lithium. 

In view of this exciting possibility, we have recently studied the equilib- 
rium properties of atomic lithium is a harmonic oscillator potential y™(x) = 
7Tia;^x^/2, with a trapping frequency of Ci;/27r — 144 Hz [Q. To incorporate 
the effect of the external potential we have, as mentioned above, applied the 
local-density approximation. The result is shown in figure 15. Note that the 
use of the local-density approximation implies that we perform the above 
outlined homogeneous calculation for each point in space with a chemical 
potential that is equal to 

1 T T ATTah n(x) 

A. X = A^ - Tjmco V Y . 2.157 

2 m 2 

The third term in the right-hand side represents the mean-field effect of the 
Hartree contribution to the selfenergy of the fermions. Since our Hubbard- 
Stratonovich procedure is in principle exact, it is not immediately clear why 
such a term must be included in the theory. It can, however, be shown that 
it arises from the fluctuations of the BCS gap parameter. With this remark, 
we end our development of the equilibrium field theory of trapped atomic 
gases. We now turn our attention to the nonequilibrium field theory. 

3 Nonequilibrium field theory 

As we have already seen several times, the equlibrium field theory gives us 
also information on dynamical properties of the gas by means of the sub- 
stitution LUn — + —iuj or equivalently t — > it. However, the Popov theory of 
Bose-Einstein condensation has shown us that we have to be very careful 
with this procedure, because we may not always end up with the correct 
physics needed for a discription of a particular experiment. It is for these 
cases that a truly nonequilibrium field theory is required and the devel- 
opment of such a theory is the topic of section p.3| . Before embarking on 
that, however, we first consider two important dynamical problems where 
the equilibrium theory does give us all the answers. 

3.1 Macroscopic quantum tunneling of a condensate 

In a Bose gas with effectively attractive interactions, i.e., with a negative 
scattering length a, a condensate will always have the tendency to collapse 
to a high density state due to the gain in energy that can be obtained in 
this way |44|. The most important question in this context is, therefore, 
if a condensate can exist sufficiently long to be experimentally observed. 
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Fig. 15. Density distribution n(x) and energy gap Ao(x) for a ^Li atomic gas 
consisting of 2.865 x 10^ atoms in each spin state a) at T = 15 nK, b) at T = 33 
nK, slightly below Tc, and c) at T = Tc = 37 nK. The left scale of each plot refers 
to the density and the right scale to the energy gap. The dotted line in c) shows 
the density distribution for an ideal Fermi gas with the same number of particles 
and at the same temperature. 



Neglecting the variation of the noncondensate density on the size of the 
condensate |^, Q, we know from the Popov theory that the dynamics of 
the collapse is, apart from an unimportant shift in the chemical potential, 
determined by the Gross-Pitaevskii equation [E5[ 



2m 



(3.1) 



+ l/°''(x)-Ai + 



Anah 



i(0(x,i))in(0(x,t)) 
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More precisely, this determines only the semiclassical dynamics. If we also 
want to study the quantum fluctuations, which is necessary if we are also 
interested in how the condensate tunnels through the macroscopic energy 
barrier, it is most convenient to calculate the grand canonical partition 
function of the condensate ^] . Quantizing the Gross-Pitaevskii equation 
we obtain for this partition function the functional integral 

Z{p)^ j dW-]d[cj,] exp|-i5[0*,0]| , (3.2) 

over the complex field 0(x, r) and with the Euclidian action 

SW,4>]= j^ drj d^ri^,r)(h— + ^— (3.3) 

+ y-(x)-M+^^|</>(x,r)p)</>(x,r). 
m J 

As always for Bose systems, the integration is only over fields that are 
periodic on the imaginary time axis. 

Although it has recently been shown by Freire and Arovas that the 
tunneling process can also be studied in terms of the complex field 0(x, r) 
Ijisll , we believe that it leads to somewhat more physical insight if we use 
instead the fields p(x, r) and 0(x, r) that correspond to the density and 
phase fluctuations of the condensate, respectively. They are introduced by 
performing the canonical variable transformation ||33|| 



0(x,r) = Vp(^e^'''''"^ 
in the functional integral for the partition function. As a result we find 

Zifi) ^ f d[p]d[e] expj-i^[p,0;/i]| , (3.4) 



with 



S[p,6;fj] — dr dx i ihp{x,T) 



96i(x,- 



dr 

i2„/ . _\ t-2 



(3.5) 



+ i^^(V0(x,r))^ + ^^(Vp(x,r))^ 

+y°'^(x)p(x, r) - ^p(x, r) + ^^^p\x, t) 

m 

Next, we notice that this action is only quadratic in the phase fiuctuations. 
The field ^(x, r) can therefore be integrated over exactly, because it only 
involves the evaluation of a gaussian integral. 



H.T.C. Stoof: Field theory for trapped atomic gases 49 



Compared to ordinary gaussian integrals there is, however, one shght 
comphcation which is associated with the fact that 0(x, r) are phase vari- 
ables. This implies that the periodicity of the original field (/!)(x, r) only 
constraints the phase field 9{x, r) to be periodic up to a multiple of 27r. 



To evaluate the grand canonical partition function in equation (3.4) we 
must therefore first integrate over all fields 9{x, r) that obey the boundary 
condition 0(x, hf3) = 9{x, 0) + 2nj and subsequently sum over all possi- 
ble integers j. Because these different boundary conditions only affect the 
zero- momentum part of 0(x, r) we first have to evaluate the sum 



X] / d[eo] cxpl-i dr No{t] 



dr 



with No{t) = J dx p{x,t) the number of condensate particles. After per- 
forming a partial integration on the integral in the exponent, we can carry 
out the path integration over 9o{t) to obtain 



.2 



g-^^WoJ^ 



dNoJT) 
dr 



As expected, the integration over the global phase of the condensate leads to 
the constraint of a constant number of condensate particles, i.e., No(t) = 
Nq. Moreover, we have ^ g^'^*^"-' = J2i^i-^o — j), which restricts the 
number of condensate particles to an integer. Putting all these results to- 
gether, we see that the integration over the zero-momentum part of p(x, r) 
is only a sum over the number of condensate particles and we have that 

Zifi) = J2ef'^'''>ZN,. (3.6) 

No 

Here we introduced the canonical partition function of the condensate, 
which is apparently equal to the functional integral 

Zno ^ J d[p]d[0] exp[-^S[p,e;0]\ (3.7) 

over all the nonzero momentum components of the density and phase fields. 
The integration over the nonzero momentum components of the phase 
field 0(x, r) is easily performed, because it now involves an ordinary gaus- 
sian integral. Introducing the Green's function for the phase fluctuations 
G(x,x';p) by 

- ((Vp) • V + pV') G(x, x'; p) = 5(x - x') , (3.8) 

m 
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we immediately obtain the desired effective action for the density field 

+ J^ drj dx (^-_^_^(Vp(x, t))2 + F-(x)p(x, r) 

2'Kah 9 . 

-P (x,t) 



(3.9) 



m 

Being an action for the density fluctuations of the condensate, 5"''^ [p] also 
describes all the collisionless modes of the condensate. This is important for 
our purposes, because the mode which becomes unstable first, determines 
precisely how the condensate collapses. Moreover, it determines the proba- 
bility with which the collapse is going to take place, both for quantum and 
thermal fluctuations, since the energy barrier is smallest in that direction of 
the configuration space. It should be noted that as long as we can neglect 
the interaction between the condensate and the thermal cloud, the action 
S*^^ [p] describes also the collective modes of a gas with positive scattering 
length. For various other theoretical approaches that have been applied 
under these conditions see, for example, Refs. P9|, pOJ, P^-p5| . The actual 
measurements of the collective mode frequencies have been performed by 



Jin et al. |22] and Mewes et al. |27| and are at sufficiently low temperatures 
indeed in good agreement with the theoretical predictions [^^. We ex- 
pect the same to be true for a gas with effectively attractive interactions 
and, therefore, the action S'°*''[p] to be a good starting point for the following 
discussion. 

To obtain the collisionless modes explicitly we consider first the case of 
an ideal Bose gas by putting a = 0. For the ideal Bose gas we expect the 
gaussian profile 

1 \3/2 



p(-;'^W)-^o(^^^^j -p(^-^j (3.10) 

to describe an exact mode of the condensate. The reason is that in the 
noninteracting case we can make a density fluctuation by taking one particle 
from the condensate and putting that in one of the excited states of the 
external potential. The corresponding density fluctuation obeys 

5p(x,i)cxe-^('"-^°)*/''x;(x)xo(x). 

For the experimentally relevant case of an isotropic harmonic oscillator pl^ 
it is more convenient to use instead of the cartesian quantum numbers n, 
the two angular momentum quantum numbers £ and m and the quantum 
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number n that counts the number of nodes in the radial wave function 
Xne{x). The density fluctuation then becomes 

<5p(x,i)«e-'(2"+^)-*Xn£(x)r,:„(i)-^ 



(7r?2)3/4 ' 

with Cnim — eooo = (2?i + i)huj the excitation energy and I = {h/mui)^''^ 
the size of the condensate wave function. Comparing this now with the 
expansion of the gaussian profile in equation ( 3.10| ) around the groundstate 



density profile, which is obtained by substituting (/(t) = / + Sq{T), we find 
that 

'5P(x,t) = -V6No^—xw{x)Yo„{x)-——^ (3.11) 

has precisely the same form as a density fluctuation in which one particle 
is taken from the condensate and put into the harmonic oscillator state 
with quantum numbers {nim) = (100). The frequency of this so-called 
'breathing' mode described by the gaussian density profile must therefore 
be equal to 2uj. 

To proof that this is indeed correct, we need to evaluate the effective 
action 5'®'^[p], and hence the Green's function G'(x, x';/?), for a gaussian 



density profile. Substituting such a profile in equation (3.8) immediately 
leads to G(x, x';p) = G{x,x';q)/p{x';q), with 

-(-^K-V + vA G(x, x'; q) = 5(x - x') . (3.12) 

ra \ q-^ J 

The latter equation can be solved, if we can solve the eigenvalue problem 

V'-^^)e(x) = Ae(x). (3.13) 

This turns out to be an easy task, because substituting 

^nlrai^) ^ ^nl.{x) y(?m(x) (3-14) 

X 

gives essentially the radial Schrodinger equation for an isotropic harmonic 
oscillator with frequency ujq — fijmq^ ^ i.e.. 

The desired eigenfunctions are therefore ^n^m(x;(7) — ipnim{^)e^ , with 
¥'r!im(x) the properly normalized harmonic oscillator states with the en- 
ergies (2n + £ + 3/2)hujq, and the corresponding eigenvalues are X„e{q) = 
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—2{2n + t}/q^. Introducing finally the 'dual' eigenfunctions £,ntm{'^\<l) = 
fnimi^)^^^ ^^' ' ^^^ Green's function G{:K,x';q) is given by 



n£7n 



G(x,x';(j) = } C„^m(x;g)— ; — --6if™(x';g) . (3.16) 



Note that prime on the summation sign indicates that the sum is over all 
quantum numbers except (nim) = (000). The latter is excluded because 
the associated eigenfunction ^ooo(x;(7) is just a constant and thus docs not 
contribute to G(x, x';p), which is defined as the Green's function for all 
phase fluctuations with nonvanishing momenta. 

Putting all these results together, we see that the dynamics of the col- 
lective variable g(r) is determined by the action 




3mNo fdqV f Zr? 3 

T^ + A^o -, ^ + 7 

dT ) \ Amq^ 4 

2 



^0 1 2 + 7™^'?" (3-17) 



that is equivalent to the action of a particle with effective mass m* = 
3mNo/2 in a potential V{q) = NoiSh^/imq"^ + 3muj^q'^/4). As expected 
from our previous remarks, this potential has a minimum for q = I and can 
be expanded near its minimum as 

Viq) ~ ^Nohco + ^m*{2iu)^{dqf . (3.18) 

It thus comfirms that the gaussian profile describes a breathing mode with 
frequency 2a; around an equilibrium density profile that is given by p(x; /) = 
A^o 1X000 (x)|2. 

Our next task is to investigate how interactions affect this result. Con- 
sidering again only gaussian density profiles, the action S'°'^[g] is again that 
of a particle with effective mass m* — 3mNo/2 but now in the potential |5^ 

V{q) = No —^ + -mu;V - ° '3 • (3.19) 

The physically most important feature of this potential is that it is un- 
bounded from below, since V{q) -^ —00 if g | 0. Hence, the conden- 
sate indeed always has the tendency to collapse to the high-density state 
limqio p(x; q) = NoS{x.). However, if the number of condensate particles is 
sufficiently small, or more precisely if p7[ 

2\/2^ II , ^ 
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the condensate has to overcome a macroscopic energy barrier before it 
can collapse. Under these conditions the condensate is therefore really 
metastable and can in principle be observed experimentally. The most im- 
portant question in this respect is of course: How metastable is the conden- 
sate? Within the gaussian approximation this question is easily answered, 
because then the dynamics of the condensate is equivalent to the dynamics 
of a particle in an unstable potential, as we have just seen. We therefore 
only need to evaluate the WKB-expression for the tunnneling rate [Q and 
compare this to the rate of decay due to thermal fluctuations by calculating 
also the height of the energy barrier. The outcome of this comparison for the 
conditions of the experiment with atomic ^Li is presented in reference IM] 
and shows that, for the relatively high temperatures T ^ huj/kB that have 
been obtained thusfar |0] , the decay by means of thermal fluctuations over 
the energy barrier is the dominant decay mechanism of the condensate. 

More important for our purposes, however, is that sufficiently close to 
the maximum number of condensate particles iVmax the collective decay 
of the condensate discussed above is always much more probable than the 
decay due to two and three-body collisions that lead to a spin-flip or the 
formation of ^Li molecules, respectively. As a result the collapse of the 
condensate should be observable within the finite lifetime of the gas. In 
fact, on the basis of this separation of time scales we expect the condensate 
to go through a number of growth and collapse cycles |£^, Q . Physically 
this picture arises as follows. Starting from a gas with a number of atoms 
N 3> A'max, the condensate will initially grow as a response to evaporative 
cooling. However, if the number of condensate atoms starts to come close 
to A^max, the condensate fluctuates over the energy barrier and collapses 
in a very short time of 0(l/w) pl[. During the collapse the condensate 
density increases rapidly and two and three-body inelastic processes quickly 
remove almost all the atoms from the condensate. After this has occurred 
the condensate grows again from the noncondensed part of the gas and a 
new growth and collapse cycle begins. It is only after many of these cycles 
that enough atoms are removed for the gas to relax to an equilibrium with 
a number of condensate particles that is less than A^max- This is shows 
quantitatively in figure 16 for the experimental conditions of interest. 

A final issue which needs to be addressed at this point is the actual 
dynamics of the collapse and, in particular, how we must include the effect 
of the inelastic growth and decay processes on this dynamics. Unfortunately, 
the inclusion of these effects is rather complicated and we need the results 



of section 3.4 to be able to correctly address this problem. It is, however, 
of some interest because Sackett et al. have recently observed that after a 
single collapse there remains a remnant of the condensate with about 10% 
of the initial number of atoms lO. At present, it is an important open 
problem to theoretically understand the magnitude of this remnant. 
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Fig. 16. Typical evolution of condensate number A'o in response to evaporative 
cooling. The long-time decay of the condensate is due to two and three-body 
inelastic collisions. 



3.2 Phase diffusion 



As we have just seen explicitly, a particularly interesting consequence of 
the finite size of the gas is that quantum fluctuations play a much more 
important role. Although this is especially true for the case of attractive 
interactions that we considered in section 3T , it is also true for a Bose gas 
with repulsive interactions. A striking example in this respect is the phe- 
nomenon of phase 'diffusion', which was recently discussed by Lewenstein 
and You |^2[. We rederive their results for a trapped Bose gas in a moment, 
but first consider also the same phenomenon for a neutral and homoge- 
neous superconductor. In this manner it is possible to bring out the physics 
involved more clearly. 



Using the approach of section 2.8.2, it can be shown that at zero tem- 
perature the dynamics of the superconducting order parameter, i.e., the 
BCS gap parameter A(x, t) that is proportional to the wave function of the 
condensate of Cooper pairs, is in a good approximation determined by a 
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time-dependent Ginzburg-Landau theory |^,^,Q with the action 

f n' 



S^-[A*,A] = ^JdtJ. 



|AoP 



dA ^ 



dt 



(3.21) 



Si.\VA\^ + 2\A\-(l-M 




where N{0) is the density of states for one spin projection at the Fermi 
energy ep = rnvp/2 and Ag is the equihbrium value of the order parameter 
pql . Writing the complex order parameter in terms of an amplitude and a 
phase, we immediately observe that the amplitude fluctuations are gapped 
p6[ | and can, therefore, be safely neglected at large length scales. The long- 
wavelength dynamics of the superconductor is thus dominated by the phase 
fluctuations, according to the action 



(3.22) 



This also implies that the global phase 9o{t) — J dK 9{x, t)/V of the super- 
conductor has a dynamics that is governed by 

using the fact that the total volume V of the system is given by N/n. 

Up to now our discussion has again been semiclassical. To consider also 
the quantum fluctuations, we have to quantize this theory by applying the 
usual rules of quantum mechanics. Doing so, we find that the wave function 
of the overall phase obeys a Schrodinger equation 



9 n 9^ 

^81^°'*'^ 'mm del 



in-^i0o; t) = -^^j;^ T^*(^o; i) , (3.24) 



with a 'diffusion' constant that can easily be shown to be equal to the 
quantity [2/h)dep/dN |67|] and is, most importantly for our purposes, pro- 
portional to 1/A^. In the thermodynamic limit A^ -^ oo a state with a well 
defined stationary phase is clearly a solution and we are then dealing with a 
system having a spontaneously broken U{1) symmetry. However, for a finite 
(and fixed) number of particles the global phase cannot be well defined at all 
times and always has to 'diffuse' in accordance with the above Schrodinger 
equation. Note also that in the groundstate the phase is fully undetermined 
and |^(6'o;t)p = l/27r. Maybe surprisingly, the same calculation is some- 
what more complicated for a Bose gas because the amplitude fluctuations of 
the order parameter cannot be neglected even at the largest length scales. 
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However, taking these amphtude fluctuations into account properly, we nev- 



ertheless arrive at an action that is equivalent to equation (3.23) and hence 
again leads to the phenomenon of phase 'diffusion'. 

We start again from the action S[p, 9; /x] for the condensate. The differ- 
ence with the previous subsection is, however, that now we are not so much 
interested in the dynamics of the density but in the phase dynamics instead. 
Therefore we now want to integrate over the density field p(x, r) . This can- 
not be done exactly and we therefore here consider only the strong-coupling 
limit, which was also treated by Lewenstein and You p2[ . In that limit we 
are allowed to neglect the gradient of the average density profile |6^ and 
the action S[p, 9] /i] is for the longest wavelengths well approximated by 

S[p,9;p.] = r dT f dx ftnp{x,T)^^^^^ (3.25) 

+ \^^'^(x)p(x, r) - MP(X, r) + ^^^p\x, t) 

m 

In equilibrium the average density profile of the condensate thus obeys 

(p(x)) = -^ {p - y-(x)) e{p - y-(x)) . (3.26) 

Anan 

Performing now the shift p(x,T) ~ (/5(x)) + (5p(x, t), we find for the zero- 
momentum part of the action p9] 

S[dNo,9o;p] = nf3Eo{p)+ f dr UhSNo^ + ^^^ {6Na)A , (3.27) 

where Eo{p) and Vo{p) correspond, respectively, to the energy and the 
volume of the condensate in the so-called Thomas- Fermi approximation [Q . 
Moreover, 6No{t) = J dx Jp(x, r) represents the fluctuations in the total 
number of condensate particles in that same approximation, implying that 
the density fluctuations (5p(x, r) are only nonzero in that region of space 
where the condensate density does not vanish. 

Performing now the integration over the number fluctuations SNo(t) 
and the usual Wick rotation to real times r — )■ it, we immediately see that 
the effective action for the global phase of the condensate has precisely the 



same form as in equation (3.23), i.e 



-cs,n . ..1 _ mVo{p) /■ j^ fd9o 



The appropriate 'diffusion' constant is therefore equal to 2T:ah/mVo{p), 
which can easily be shown to be equal to {l/2fi)dp/dNo if we make use 
of the fact that in the Thomas- Fermi approximation the chemical potential 
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obeys fi = muj'^R'^Tp/2 and the radius of the condensate is given by Rtf = 
{15ah? No/m'^uj'^)^''^ Q. Hence, the 'diffusion' constant is proportional to 

I/NJ . Note that if the condensate where contained in a box the 'diffusion' 
constant would be proportional to 1/Nq instead, ft is important to note 
also that, in contrast to the case of a fermionic superfluid, we have to 
integrate over the amplitude fluctuations of the order parameter to arrive 
at a quadratic action for the phase fluctuations. This leads to the important 
conclusion that for a bosonic superfluid it is impossible to be in a state with 
only phase fluctuations and no density fluctuations, even at the largest 
length scales. 

3.3 Quantum kinetic theory 

For simplicity we restrict ourselves for the rest of this course to the case 
of a doubly-polarized atomic Bose gas, which for most practical purposes 
can be seen as a gas of spin-less bosons Q. Once we have understood 
the quantum kinetic theory for this particular case, however, we can easily 
generalize to multi-component Bose gases that are receiving considerable 
attention at present ItI]. Furthermore, the same methods can also be used 
for fermionic gases and have recently already been applied to the problem 



of Cooper-pair formation in a two-component fermion gas 1 73 



3.3.f Ideal Bose gas 

In textbooks an ideal Bose gas is generally discussed in terms of the average 
occupation numbers of the one-particle states Xn (x) fTJ] . Given the density 
matrix p(io) of the gas at an initial time t^, these occupation numbers obey 

iV„(t) = Tr [p(io)^i(i)^„(t)] , (3.29) 

with 'ipl^{t) and V'n(i) the usual (Heisenberg picture) creation and annihila- 
tion operators of second quantization, respectively. Because the hamiltonian 
of the gas 

H = Y,e^{^i{t)Mt) (3.30) 

n 

commutes with the number operators iV„(t) = ^^(t)^n(i), the nonequilibri- 
um dynamics of the system is trivial and the average occupation numbers 
are at all times equal to their value at the initial time t^. If we are also 
interested in fluctuations, it is convenient to introduce the eigenstates of the 
number operators, i.e.. 
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and to consider the fuU probabiUty distribution 

Pi{N^};t)^Tv[p{to)\{N^};t){{N^};t\] , (3.32) 

which is again independent of time for an ideal Bose gas. The average 
occupation numbers are then determined by 

Nn{t) = Y. A^„F({iVn};i) (3.33) 

{JV„} 

and the fluctuations can be obtained from similar expressions. 

As indicated in section |3.2| , we are also interested in the phase dynamics 
of the gas. In analogy with the occupation number representation in equa- 



tion (3.32), we can now obtain a completely different description of the Bose 
gas by making use of coherent states and considering the probability distri- 
bution P[(f)* ,(j)]t\. Although we expect that this probability distribution is 
again independent of time, let us nevertheless proceed to derive its equa- 
tion of motion in a way that can be generalized lateron when we consider 
an interacting Bose gas. First, we use that 

Pr,0;t] = f d[r,]d[M P[ro,^oM ,^,l[^]^]i"'l']^' ,. ■ (3.34) 
J {(P,tmt){<po;to\(po..to) 

This is a particularly useful result, because the time dependence is now 
completely determined by the matrix element ((/); t|(/)o;^o) for which the 
functional integral representation is well known from our discussion in sec- 
tion |2.3| . It is given by a 'path' integral over all complex field evolutions 
V'(x, f+) = J2n V'n(i+)Xn(x) from to to t. More precisely we have 

V*(x,t)=0*(x) c ■ 



(0;i|0o;to)= / d[r]dm e^p\-S+[r,i^]> , (3.35) 

with the forward action S+[ip*,'iJj] given by 

S+[r,i'] = E {~^Wnmn{t) (3.36) 

+ ^ dt+ rn{t+) [ih^ - en) Mt+)\ ■ 

In the same manner the matrix element {(f>;t\(j)Q;tQ)* = {(/jq', tolcf); t) can be 
written as a 'path' integral over all possible field configurations ip{'x.,t^) — 
X^n V'n(i-)Xn(x) cvolving backward in time from t to t^, i.e., 

<.V*(x,to)=</'S(x) r : > 

(0;#o;io>*=/ d[r]dm exp\-S.[r,i']\ , (3.37) 

J i>ix.t)=<p{x) L " J 
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with a backward action 



5_[^*,V]=5]|-*?iC(io)V'n(io) 



(3.38) 



+ / di_ ^;(t_) U?i— -enj^n(i_) 

n ^ 

+ l" dt^ ^„(t_) (-ih^ - en) rn{t-)\ ■ 

Putting all these results together, we see that the probability distribution 
P [(/)*, cj);t] can in fact be represented by a functional integral over all fields 
ip{x,t) that evolve backwards from t to to and then forward in time from to 
to t. Absorbing the factor P[(l)Q,(j>o;to] into the measure of the functional 
integral, we thus arrive at the desired result that 

7A*(x,t)=</.*(x) f ■ 



PWA-A= diV'IdM exp -5[V;*,V] , (3.39) 

JV(x,t)=0(x) I "- J 

where the total (backward-forward) action in first instance obeys 

s[r,^] = s-[r,^ + s+[r,^] = -inY, (rnWMt) - i^i^nf) (3.40) 

n 

- enC(t')^n(i') 



and the integration along the Schwinger-Keldysh contour C* is defined by 
f^t dt' = J*° dt_ + J^^ dt+ |f|-0. Note also that in equation ( |3.39D we 



have explicitly specified the boundary conditions on the functional integral. 
It is interesting to mention that these boundary conditions are essentially 
dictated by the topological terms in the action 5'[V'*, V']j which is a general 
feature of the path-integral formulation of quantum mechanics due to the 
fact that the quantum theory should have the correct (semi)classical limit 
jTSll . Making use of the periodicity of the field '0(x, t) on the Schwinger- 
Keldysh contour, the variational principle ^5'['0*,'0]/^'0n(^±) — indeed 
leads not only to the Euler-Lagrange equation 

in-^Mt±) = enV'n(t±) , (3.41) 
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which agrees with the Heisenberg equation of motion i?idipn(t±)/dt± = 
[il'n(t±),H] for the annihilation operators, but also with the appropriate 
boundary condition 6^p^{t) — 0. In the same way we find the complex 
conjugate results if we require th at SS lip* ,ip]/S'ipn{t±) = 0. Substituting the 
boundary conditions in equation (3.40) and performing a partial integration, 
we then finally obtain for the action 



5[V'*,V]=5]£di'C(0(*?i|7 



^n(t') , 



which we for completeness sake also rewrite as 






(3.42) 



(3.43) 



2m 



-y"''(x)U'(x,i') 



We are now in a position to derive the equation of motion, i.e., the 
Fokker-Planck equation, for the probability distibution P[(j)* ,(f>;t]. This is 
most easily achieved by performing the variable transformation ■(/'(x, t±) = 
(j){x, t')±£^{K, t')/2 in equation ( 3.3£ ). In this manner the fields ^/'(x,i_) and 
■(/'(x, i_|_) that live on the backward and forward branch of the Schwinger- 
Keldysh contour, respectively, are 'projected' onto the real time axis. More- 
over, at the same time we effect a separation between the (semi)classical 
dynamics described by (/)(x, i') and the quantum fluctuations determined 
by f (x, t'). After the transformation we have 



r.0'(x,t)=0*(x) 

I 

;x,t)=0(x) 



(3.44) 



w]d[^]ld[ew] 



exp<j-5[^*,0;r,C] 



with 



^r,'/';r,C] = E/ ^*' {C(i')(»?i^-en)cn(i') (3.45) 

+ C(O(*?i|7-en)0n(O} • 

Because this action is linear in ^n(i') and £,^{t'), the integration over the 
quantum fluctuations leads just to a constraint and we flnd that 



PI 



■,t\ = 



0*(x,t)=0'(x) 
<x,t)=0(x) 



(3.46) 



-4-TVn(^')-f4-T)^"(0 



dt' 



dt' 
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or equivalently that 

P[0*,0;t]= /'d[05]d[0o]P[0S,'^o;io]nMl'^n-</'^„'(i)l') , (3.47) 



where we introduced the quantity cj)'^ (t) obeying the semiclassical equation 
of motion 

z?i^(/.J,\i) = 6„C'W (3.48) 

and the initial condition 0n(*o) — 0O;n. 

The latter equation is thus solved by (f>n{t) — (/)o;ne^"^"^*^*''-'/'' and we 



conclude from a simple change of variables in equation (3.47) that for an 
ideal Bose gas P[(j)*,(j>;t] — P[0*,(/);io], as expected. We also see from 
equation ([3.47) that the desired equation of motion for P[(j)* , (j); t] reads 



in^Plr,<f>;t] = - I y^^t„4,„ I p\4.',4,;t] (3.49) 






[(/)*, 0;t] 



This is indeed the correct Fokker-Planck equation for an ideal Bose gas [ra . 
To see that we first consider ((/'n)(0 — J d,[(l)*]d[(j>] 0„F[(/)*, 0; t]. Multiplying 
equation (|]49) with </)„ and integrating over 0(x), we easily find after a 
partial integration that 

in^A4'^){t)^en{K){t) , (3.50) 

ot 

which precisely corresponds to the equation of motion of the expectation 
value (■0n(i)) = Tr[p(io)V'n(0] ^^ the operator formalism. Similarly, we find 
that 

z?.|(C)(i) = -^n(C>(0, (3.51) 

in agreement with the result for (V^i(i)) = Tr[p(io)'0i(i)]. 

Next we consider the average of l^nP, for which we immediately obtain 

ih^^ilKl'm ^ 0. (3.52) 

We expect this result to be related to the fact that in the operator formalism 
the occupation numbers Nn{t) are independent of time. Although this turns 
out to be true, to give the precise relation between (|</)nP)(i) and Nn{t) 
is complicated by the fact that at equal times the operators ipnit) and 
tpiiit) do not commute. However, the path-integral formulation of quantum 
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mechanics is only capable of calculating time-ordered operator products [Q 
In our case this implies that (|0nP)(i) is the value at t' = t of 



Tr 



p{to)Ta [i^^mlit'))] (3.53) 

= e(t,i')Tr [p(to)^„(t)^i(t')l + e(i',t)Tr [p(io)^i(i')^n(i)l , 



with Tgt the time-ordening operator on the Schwinger-Keldysh contour and 
Q{t,t') the corresponding Heaviside function. Since the Heaviside function 
is equal to 1/2 at equal times pl|, we conclude that 

(|0nP)(i)=iVnW + ^ (3.54) 

and that equation ( [3.52| ) is thus fully consistent with the operator formalism. 
An intuitive understanding of the relation between (|0nP)(^) and the occu- 
pation numbers can be obtained by noting that, since all our manipulations 
up to now have been exact, we expect (|(/'nP)(i) to contain both classical 
and quantum fluctuations. These correspond precisely to the contributions 
Nri{t) and 1/2, respectively. 

Finally we need to discuss the stationary solutions of the Fokker-Planck 
equation. It is not diflicult to show that any functional that only depends on 
the amplitudes \(i>n\'^ is a solution. As it stands the Fokker-Planck equation, 
therefore, does not lead to a unique equilibrium distribution. This is not 
surprising, because for an isolated, ideal Bose gas there is no mechanism for 
redistributing the particles over the various energy levels and thus for relax- 
ation towards equilibrium. However, the situation changes when we allow 
the bosons in the trap to tunnel back and forth to a reservoir at a tempera- 
ture T. The corrections to the Fokker-Planck equation that are required to 
describe the physics in this case are considered next. However, to determine 
these corrections in the most convenient way, we have to slightly generalize 
the above theory because with the probability distibution P[4>* , 4>] t] we are 
only able to study spatial, but not temporal correlations in the Bose gas. 

To study also those we follow the by now well-known procedure in quan- 
tum field theory and construct a generating functional Z[J, J*] for all (time- 
ordered) correlation functions. It is obtained by performing two steps. First, 
we introduce the probability distribution Pj [cj)* , (/); t] for a Bose gas in the 
presence of the external currents J(x, t) and J*(x, t) by adding to the hamil- 
tonian the terms 

h / dx ('^(x,i)J*(x,i)-F J(x,i)?/;t(x,i)') = 

-n^(^„(i)J:(i) + Jn(i)^i(0 
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As a result we have 
Pj[^\^-t]^ d[r]dme^I>{TS[r,^]\ (3.55) 

"'i/'(x,t)=0(x) L" J 

xexpii / dt' I dx {tP{x,t')J*{x,t') + J{x,t')tP*{x,t')) 

Second, we integrate this expression over (/)(x) to obtain the desired gener- 
ating functional. Hence 

Z[J,J*] = j dW]d[^] PjW,M = fd[r]dm expUs[r,i']\ (3.56) 



xexpii dt dx {ipix, t)J* (x, t) + J(x, t)^* (x, t)) 

It is important to realize that Z[J, J*] is indeed independent of the time 
t because of the fact that Pj [0* , (f)\ t] is a probability distribution and thus 
properly normalized. We are therefore allowed to deform the contour C* to 
any closed contour that runs through t^. Since we are in principle inter- 
ested in all times t > to, the most convenient choice is the countour that 
runs backward from infinity to to and then forwards from to to infinity. 
This contour is denoted by C°° and also called the Schwinger-Keldysh con- 
tour in the following because there is in practice never confusion with the 
more restricted contour C* that is required when we consider a probability 
distribution. With this choice it is also clear that Z[J,J*] is a generating 
functional. Indeed, equation ( p.56|) shows explicitly that all time-ordered 
correlation functions can be obtained by functional differentiation with re- 
spect to the currents J{x,t) and J*{x,t). We have, for instance, that 



Tmom^,t)]^]jj^^zij,r] 



(3.57) 

J,J*=0 



and similarly that 



Tr 



■p(io)rc^(^(x,t)V3t(x',i'))] (3.58) 

Z\J,J* 



i^6J*ix,t)6Jix',t') 



J,J'=0 



Note that the times t and t' always have to be larger or equal to to for these 
identities to be valid. 

3.3.2 Ideal Bose gas in contact with a reservoir 

Inspired by the Caldeira-Leggett model for a particle experiencing friction 
P4I , we take for the reservoir an ideal gas of N bosons in a box with volume 



64 Les Houches summer schools 



V. The states in this box are labeled by the momentum Kk and equal to 
Xk(x) = e^^'^/^/V. They have an energy e(k) = ?i^k^/2m + A^'', where 
AV'^'^ accounts for a possible bias between the potential energies of a particle 
in the center of the trap and a particle in the reservoir. The reservoir is also 
taken to be sufficiently large that it can be treated in the thermodynamic 
limit and is in an equilibrium with temperature T and chemical potential fj, 
for times t < Iq. At ip it is brought into contact with the trap by means of 
a tunnel hamiltonian 

^'"' = ^EE {tnmnml{t)+tuk)MmUt)) , (3.59) 

with complex tunneling matrix elements in(k) that for simplicity are as- 
sumed to be almost constant for momenta hk smaller that a cutoff hkc but 
to vanish rapidly for momenta larger than this ultraviolet cutoff. Moreover, 
we consider here only the low-temperature regime in which the thermal de 
Broglie wavelength A = {27rh /mkBTY'^ of the particles obeys A ^ 1/kc, 
since this is the most appropriate limit for realistic atomic gases. 

To study the evolution of the combined system for times t > t^ we thus 
have to deal with the action 

S[r,t, rn, ^b] ^Y.Jj^ ^n(*) (^^^ - en + m) V^n(t) (3.60) 

^T.j^J^ V^k(0 (*?i^ - e(k) + m) Mt) , 

-^T.T.1 ^^ (t„(k)Vn(i)V'kW+c(k)V'kWcw) 

if we measure all energies relative to the chemical potential and also intro- 
duce the complex field i/'fl(x, t) = '^^ i/)k(Oxk(x) for the degrees of freedom 
of the reservoir. However, we are only interested in the evolution of the Bose 
gas in the trap and therefore only in the time-ordered correlation functions 
of this part of the system. The corresponding generating functional 

Z[J,r] = j d[r]d[il>] j d[rMM expj j5[r,^;V'fl,^fl]| (3.61) 

xexpii / dt dx (■0(x, t) J*(x, t) -f J(x, t)V'*(x, t)) 



is of the same form as the functional integral in equation ( 3.5(]| ) , but now 
with an effective action that is defined by 

expUs°«[r,i^]\ = f d[rR]d[i^R] expUs[r,i';rR,i'R]] ■ (3.62) 
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Hence, our next task is to integrate out the field V'fl(x, t), which can be 
done exactly because it only requires the integration of a gaussian. 

To familiarize ourselves with the Schwinger-Keldysh formalism, we per- 
form the gaussian integration here in some detail. In principle, this can of 
course be done explicitly by making use of the fact that the initial density 
matrix pR{tQ) of the reservoir can be expanded in the projection operators 
on the coherent states by using the coefficients 

Pi.[|<^«P;to]=^]^^~''^''/"'^'^ (3-63) 

with A^(k) = l/(e'^(^(^)^^) — 1) the appropriate Bose distribution function 
and (3 = l/ksT. However, in practice it is much more easy to use a different 
procedure. It is based on the observation that if we introduce the S function 
on the Schwinger-Keldysh contour defined by L^ dt' 6{t, t') = 1 and the 
Green's function G(k;t,t') obeying 

ih^ - e(k) + Ml G(k; t, t') = hSU, t') , (3.64) 

at J 

the action S'[V'*, V'! V'/j? V'k] can be written as a complete square, or more 
precisely as the sum of two squares 5*1 [i/'*, "0] a-nd S'2['!/'*, V"; '^*ri V'fl:] that are 
given by 

SiWM = Y.jj^ ^n(i) U^^^ - £„ + A^) V'n(t) (3.65) 

-^J2J2[ ^*/ di'V;(i)C(k)G(k;i,t')in'(k)V'n'(i') 

and 

S2[r .^-.rR.^R] (3.66) 

= J2 [ '^* ( ^k(t) - r^fE / ^*' C(k)C(i')G(k;t',t) I 
X («?i^-e(k)+M 

respectively. Since the first term is independent of the field ipfi{'x.,t), we 
only need to evaluate J d['ip'^]d['ipfi] cxp (iS'2[V'*, V'! V'lj: V'-r]/^)- Performing 
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a shift in the integration variables, we however see that this functional 
integral is equal to 

dirnH^R] exp U ^ y ^ dt V-i^lt) U^^ - e(k) + m) Mt) \ (3.67) 

= TrK(to)] = 1 . 

As a result the effective action S"^ [ip* , ifi] is just equal to Si[ip*,ip], which 
can be slightly rewritten to read 

S'«ff[^*>] = ^ /" dt [ dt' (3.68) 



n,n' *^ ^ 



X V;(i) I Un^^ -en+fA Sn,n'Sit,t') - ;iS„,„, (t, t') | V^n'(i') , 

with a selfenergy I]„^n'(*i*') obeying 

?iS„,„'(i>i') = ^E*nW^(^'*'*')*»'W • (3.69) 

k 

This is our first example of an effective action describing the nonequi- 
librium dynamics of a Bose gas. Before we can study its consequences we 
clearly first need to determine the Green's function G(k;t,t') in terms of 
which the selfenergy is expressed. Although we know that this Green's func- 



tion fulfills equation (3.64), we cannot directly solve this equation because 
we do not know the appropriate boundary condition aX t = t' . To calculate 
G(k;t,t') we therefore have to proceed differently. It is however clear from 



equation (3.64) that G(k; t, t') is a property of the reservoir and we thus ex- 
pect that it can somehow be related to a time-ordered correlation function 
of this rcscrvior. To see that explicitly we consider again the generating 
functional of these correlation functions, i.e., 

ZR[J,.r] (3.70) 

= j d[rR]d[M exp U E /^ d^ ^k(^) (*^^ - <'') + 1^ v^k(i) I 

xexpiz^y^dt {Mi)JL{t) + MtW^{t))\ ■ 

It is again a gaussian integral and can thus be evaluated in the same manner 
as before. The result is now 



Zj7[J,J*] =exp<^-iy / dt dt' J^{t)G{k;t,t')Jk{t')} , (3.71) 
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which shows by means of equation (3.5S) that 



iG(k;i,t') = Tr [pR{to)Tc^ (^kWV^k(i') 



(3.72) 



Note first of all that this is indeed consistent with equation (3.64), be- 
cause the right-hand side obeys 



'^Wt^' 



PR{h)Tc^ (^k(i)V^I(i') 
iM(t,f')Tr [pfl,(io)[^k(t)>I(i)] 
-HTr 



(3.73) 



ihd{t,t') + (e(k) - fi)TT U(io)Tc=o (^k(i)^i(i': 



in the operator formalism. Moreover, from this identification we see that the 
desired solution fulfilling the appropriate boundary conditions is apparently 

G{k;t,t') (3.74) 

^ _jg-»(E(k)-AO(t-t')A {e(i^ i') (1 + 7v(k)) + e(i', t)N{k)} . 

The specific dependence on the backward and forward branches of the 
Schwinger-Keldysh contour is thus solely determined by the Heaviside func- 
tion 8(t,i'). As a result it is convenient to decompose the Green's function 
into its analytic pieces G^ (k; t — t') and G^ (k; t — t') [p3| by means of 



G(k; t, t') = e(i, t')G> (k; t^t') + 0(1', <)G<(k; t - t') 



(3.75) 



Due to the fact that we are always dealing with time-ordered correlation 
functions, such a decomposition turns out to be a generic feature of all the 
functions on the Schwinger-Keldysh contour that we will encounter in the 
following ^J]. For a general function F{t,t'), it is however also possible to 
have (5-function singularities. If that happens the correct decomposition is 

F{t,t') = F\t)6{t,t') + e{t,t')F> {t,t') + Q{t' ,t)F<{t,t') . (3.76) 

This more general decomposition is not required here, but will be needed in 



3.4 when we determine the effective interaction between two atoms 



section 
in a gas. 

Having obtained the Green's function of the reservoir, we can now re- 
turn to our discussion of the effective action S'^^[tp*,^p] for the Bose gas 
in the trap. Although we have chosen to derive in equation (3.68) the ef- 
fective action for the generating fmictional Z[J,J*], it is straightforward 
to show that the effective action for the probabihty distribution P[(j)*,(p; t] 
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is obtained by only replacing the contour C°° by C*. In the following we 
therefore no longer always specify all the boundary conditions on the time 
integration, if the precise details of this integration are not important and 
the discussion applies equally well to both cases. Keeping this in mind, we 
now again perform the transformation ■0n(^±) = 0n(^) i ^„(t)/2 to explic- 
itly separate the (semi)classical dynamics from the effect of fluctuations. It 
leads in first instance to |g^ 

^°*^[0*,0;r,e]-E/ dtKit)(iri^^-^n + fi]ut) (3.77) 

+ E / ^* ^n(i) (i^^i - en + A Mt) 

-J2 f dt f dt' Km^i:Ut-t')^r.'it') 

nW-Jto J to 

-J2 [ dt f dt' cm^iUt~t')^r.'it') 

where we introduced the retarded and advanced components of the sclfcn- 
ergy 

^tlit - t') = ±Q{±{t - t')) (E>„,(t - t') - S]<„,(t - t')) (3.78) 

that affect the terms in the action that are linear in ^*(t) and i,■n.{t)^ respec- 
tively, and also the Keldysh component 

^lAt t') = (E>„,(t - t') + I]<„,(i - t')) (3.79) 

that is associated with the part quadratic in the fluctuations. 

The physical content of these various components of the selfenergy is 
understood most clearly if we now apply again the beautiful functional- 
integral procedure that is due to Hubbard and Stratonovich. The basic 
idea is to write the factor 



^^p^-^E 



Y. I dt I dt' cAmlAt^t')s,^'{t')\ 

n,n'-^*o Jto J 



in the integrant of / d[0*]d[^] / d[£;]d[i] exp {iS^^lcj)* ,(j); T, 'C]/^) as a gaus- 
sian integral over a complex field ?7(x, i). It is equivalent, but in practice 
more convenient, to just multiply the integrant by a factor 1 that is written 
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as a gaussian integral J d[ri*]d[r]] exp (iS'®'^[?7*, TyJ/Ti) with 

S'^'iv^v] (3.80) 

= ^J2 [ dt f dt' [ 2ry;(i) -J2 f dt" C"inJ^^n",nit" - t)] 

a complete square. Adding this to S°^[(j)* , (p; £,*,£] the total effective action 
becomes 

5^*^[0*,(/.;r,^;^*,^] (3.81) 

= Y,jdtjdt' cj,l{t) I Un^^ - e„ + A* - 77; (i)") (5„,„,^(i - t') 

-?.<2,(t-t')}en'(t') 

+ ^ I dt I dt' Cit) I h?J^ - Cn + /^ - 77n(i)j ^n,n"5(i " t') 

-?iEW(i-t')|0„'(O 
+2^ /■ d^ /■ dt' r7;(i)(?iS^')„i,(t-t')r;„,(i') 

and is thus linear in ^n(0 ^^'^ Cn(0- Integrating over these fluctuations 
we conclude from this action that the field 0(x, t) is constraint to obey the 
Langevin equations 86 

o poo 

*?i7^0n(i) = (e„-Ai)0n(t) + V / dt' h^'^+l,{t-t')(br,it')+7^^{t) (3.82) 
and 

o pOG 

-in-6l{t) = (e„-^)0;(t)+^ / dt' <f>l,{t')h^i';lit'-t)+r,Ut) (3.83) 
with gaussian noise terms ?7n(i) and ??*(i) that from the last term in the 



right-hand side of equation (3.81) are seen to have the time correlations 



{v:(.t)vn'it')) = '-^^^,At-t') (3.84) 

= lIWf (l + 2iV(k))C(k)e-(^('')-'')(*-*')/^'i„,(k) 
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in the thermodynamic limit. Moreover, in that limit the retarded selfenergy 
is equal to 

?^S1^^'(^-^') =-e(t-t')j/ (0J C(k)e-'(^('')"'')(*-*')/''t„,(k) (3.85) 

and the advanced selfenergy can be found from the relation E^, [^(t' — t) = 
( E^ j^, (t — t') J . Note that the Heavisidc functions in these retarded and 



advanced selfenergies are precisely such that Eqs. (3.82) and (3.83) are 
causal, i.e., the evolution of 0(x, t) and (/>*(x, t) depends only on their value 
at previous times. This is clearly a sensible result and explains why these 
components of the full selfenergy enter into the terms in the action that are 
linear in the fluctuations. 

To understand why the Keldysh component enters into the terms that 
are quadratic in the fluctuations is more complicated and is related to the fa- 
mous fluctuation-dissipation theorem p3| , |8q , |8^ . Physically, the fluctuation- 
dissipation theorem ensures in our case that, due to the interactions with 
the reservoir, the gas in the trap now relaxes towards thermal equilibrium in 
the limit t —^ oo. Let us therefore study this limit in somewhat more detail. 
However, before we do so we would like to mention for completeness that if 
we apply the above Hubbard-Stratonovich transformation to the functional 
integral for the probability distribution P[(f>* , (j)] t],we immediately find that 



analogous to equation (3.47) 



P[0*,0;t] = / d[r^*]d[fj] P[fj\r^] / d[0*]d[0o] P[|0o|';io] (3.86) 



where (/)„(i) now obeys the Langevin equation in equation (3.82) for a par- 
ticular realization of the noise rin{t) and with the initial condition (f)'^{to) — 
0O:n- Furthermore, the probability distribution for each realization of the 
noise is given by 

(3.87) 

to J to 

This result shows explicitly that, as expected, the Fokker-Planck equa- 
tion associated with the Langevin equations in Eqs. ( ^.82 ) and ( 3.83| ) is 




in fact the desired Fokker-Planck equation for the probability distribution 
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After this shght digression, we now return to the long-time behaviour of 
the gas. In the long-time limit two important simplifications occur. Looking 
at the expresions for the selfenergies, we see that their width as a function of 
the time difference t—t' is at most of order 0{Ti/kBT) in the low-temperature 
regime of interest where the thermal de Broglie wavelength of the particles 
obeys A :^ 1/fcc- If therefore t ^ to + 0{h/kBT) we are allowed to take 
the limit to -^ — oo. Taking this limit means physically that we neglect 
the initial transients that are due to the precise way in which the contact 
between the trap and the reservoir is made, and focus on the 'universal' 
dynamics that is independent of these details. In addition, at long times 
the dynamics of the gas is expected to be sufficiently slow that we can 
neglect the memory effects altogether. Finally, we consider here first the 
case of a reservoir that is so weakly coupled to the gas in the trap that we 
can treat the coupling with second order perturbation theory. As a result 
we can also neglect the nondiagonal elements of the selfenergies. In total 
we are then allowed to put (but see below) 

SS^^'C^ - 1') - sS^^'len - M)'5n,n'<5(t - 1') = s(,±)^^<5„,„,5(t - 1') , (3.88) 
where the Fourier transform of the selfenergies is defined by 

SS^'C* -t')-J de SS^(^) ^^^ . (3.89) 



Note that this implies in general that S^, [^ (e) — ( S^ j^, (e) ) , and in par- 
ticular therefore that Sn = ( Sn 



With these simplifications our Langevin equations become 

*?i^0n(i) = (e„ + fiS(,+) - M)'^n(t) + J?n(t) (3.90) 

and just the complex conjugate equation for cii!i(t). The retarded selfenergy 
in this equation is given by 

with e+ the usual notation for the limiting procedure £„ + *0. It clearly 
has real and imaginary parts that we denote by Sn and — i?n, respectively. 
Denoting the Cauchy principle value part of an integral by V, they obey 

^n^/(^C(k)-^t„(k) (3.92) 



and 



/dk 
^^<5(6„-e(k))|in(k)|2. (3.93) 
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The interpretation of these results is quite obvious if we consider the average 
of the Langevin equation, i.e., 

i?i^ (0„> (t) = (e„ + 5n - ii?„ - m) (M it) , (3.94) 

which is solved by 

(0„)(t) = (0„)(O)e-'(^"+^"-^)*/''e-^"*/'* . (3.95) 

Hence the real part of the retarded selfenergy Sn represents the shift in 
the energ y of s tate Xn(x) due to the coupling with the reservoir. Indeed, 
equation ( ^3.92 ) agrees precisely with the value of this shift in second order 



perturbation theory. Moreover, the fact that |(0„)(t)P = |(<^n)(0)pe~^-^"*/'* 
shows that the average rate of decay Fn of the state Xn(x) is equal to 2Rn/fi, 



which, together with equation (3.93), exactly reproduces Fermi's Golden 
Rule. 

From the equation of motion for the average {4>ii){t) we can also con- 
clude that the right-hand side of the Fokker-Planck equation contains the 
'streaming' terms 









However, there must now also be a 'diffusion' term due to the fact that the 
average {{(pnl"^) (t) is no longer independent of time since the interactions 
with the reservoir can lead to changes in the occupation numbers Nn{t). To 
obtain the 'diffusion' term we thus need to determine i'hd{\<j)n\'^){t)/dt. To 
do so we first formally solve the Langevin equation by 

0„(O (3.96) 



e 



Multiplying this with the complex conjugate expression and taking the av- 
erage, we obtain first of all 

(|0nP)(i) = e-2^"*/'' |(|</-„n(0) + '^^^Idt' e^^"*'/''} , (3.97) 



which shows that 



*?i^(|</>nP)(0 = -2*i?n(|</>nP)(i) " ^^iS^ ■ (3.98) 



H.T.C. Stoof: Field theory for trapped atomic gases 73 

Moreover, the Keldysh component of the selfenergy is given by 

/r/k 
j^{l + 27V(k))<5(6„ - 6(k))|i„(k)p (3.99) 

and therefore obeys 

aE^ = -2i{l + 2Nn)Rn (3.100) 

with JV„ = (e'^('^"~^) — 1)"^ the Bose distribution function. This is in 
fact the fluctuation-dissipation theorem, because it relates the strength of 
the fluctuations determined by hT,^ , to the amount of dissipation that is 
given by i?n- As mentioned previously, the fluctuation-dissipation theorem 
ensures that the g as re laxes to thermal equilibrium. This we can already 



see from equation ( 3.98 ), because substitution of the fluctuation-dissipation 
theorem leads to 

in^A\<l>n\^){t) = ~2iRn{\cf>n\^){t)+iRn{l + 2Nn) , (3.101) 

ot 

which tells us that in equilibrium (l^nP) = Nn + 1/2. In fact, we have 
argued that in gen eral (|0nP)(O = ^n{t) + 1/2- Substituting this identity 
in equation ( [3.101 ) and realizing that T^N^ = rniV(k) due to the energy- 



conserving S function in Fn, we indeed obtain the correct rate equation for 
the average occupation numbers 

^N^it) = -r„7V„(t) + F„iV„ (3.102) 

at 

= -rniv„(t)(i + iv(k)) + F„(i + iv„(t))iv(k) , 

that might justly be called the quantum Boltzmann equation for the gas, 
because the right-hand side contains precisely the rates for scattering into 
and out of the reservoir. 



Furthermore, equation ( 3.98 ) shows that the 'diffusion' term in the 
Fokker-Planck equation is 



because the first term in the right-hand side of equation ( [3.98 ) is also present 



without the noise and is therefore already accounted for in the 'streaming' 
terms. In total we have thus arrived at 

^n^P[r,c|,; t] = -ly] -^(e„ + SS(,+) - /i)0„ ) P[r , </-; t] (3.103) 



dt 







2 ^ d^id<t> 
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Using again the fluctuation-dissipation theorem, it is not difflcuh to show 
that the stationary solution of this Fokker-Planck equation is 

P\(j)*,(j);oo]=Y\ 5-^expj ^-^\(t>J^\ ■ (3.104) 

Although this appears to be quite reasonable, it is important to note that 
this result is not fully consistent because the 'streaming' terms in equa- 



tion (3.103) show that the energies of the states Xn(x) are shifted and equal 
to En + 5*11 ■ We therefore expect that the equilibrium occupation numbers are 
equal to the Bose distribution evaluated at those energies and not at the un- 
shifted values en. The reason for this inconsistency is, of course, that we have 
taken Sn equal to I]n,n (cn — A*), which was only justified on the basis 
of second-order perturbation theory. We now, however, see that to obtain a 
fully consistent theory we must take instead En ' = Sn.n (cn + ^n — A*) • 
This implies that to obtain the retarded selfenergy, we first have to solve 
for the real part 



dk 

(2^3 ^»--e„ + ^„_e(k)- 



Sn - / T^^TW C(k)— ^ TTut^O^) (3-105) 



and then substitute this into 

Rn = 7rJ -^ S{en + Sn- e(k))|i„(k)|2 (3.106) 

to obtain the imaginary part. The fiuctuation-dissipation theorem then 
again reads hY.^ — ~-2i(l + 2Nn)Rn but now with the equilibrium distribu- 
tion Na = (^e^(.'^"+Sn-iJ-) _ l)-i. As a consequence the quantum Boltzmann 
equation and the stationary state of the Fokker-Planck equation will now 
also contain these correct equilibrium occupation numbers [g8| . 

Having arrived at the Fokker-Planck equation of our Caldeira-Leggett- 
like toy model, we have now demonstrated the use of essentially all the 
theoretical tools that are required for a discussion of a weakly interacting 
Bose gas. Before we start with this discussion, however, we want to make a 
final remark about the effect of the nondiagonal elements of the selfenergies. 
Physically, including these nondiagonal elements accounts for the change in 
the wave functions Xn(x) due to the interaction with the reservoir. It can 
clearly be neglected if the coupling with the reservoir is sufficiently weak, 
or more precisely if |?iS^, „(en + S^ — ij.)\ is much smaller than the energy 
splitting |en' — Cnl between the states in the trap. For a quantum dot in solid- 
state phyisics |^ this is often the case and for that reason we might call our 
toy model a bosonic quantum dot. However, in magnetically trapped atomic 
gases the average interaction energy between the atoms is already slightly 
below the Bose-Einstein condensation critical temperature larger than the 
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typical energy splitting and the nondiagonal elements of the selfenergies 
are very important. Such a strong-coupling situation can, of course, also 
be treated in our Caldeira-Leggett like model. The main difference is that 
we need to expand our various fields not in the eigenstates of the trapping 
potential Xn(x) but in the eigenstates Xn(^) ^^ ^^'^ nonlocal Schrodinger 
equation 

4x:.(x) - (-^ + y'^'i^)) xUx) (3.107) 

+ J dx' Re [?iE(+) (x, x'; e'„ - /i)] x'„(x') 

where e^ are the new eigenvalues and the retarded selfenergie in coordinate 
space is ?iE(+)(x, x'; e) — X^n n' Xn(x)?iS„ „/(e)x*/(x'). In this new basis 
the nondiagonal elements of the selfenergies can now be neglected and we 
find essentially the same results as before. We only need to replace £„ + 5*11 
by e'„ and tn(k) by J2n' {I d^ Xni^lx^'i^)) in'(k). Neglecting the nondi- 
agonal elements in this basis only requires that the real part of the retarded 
selfenergy is much larger that its imaginary part, which is always fulfilled 
in our case because fee A ^ 1. 

Summarizing, the coherent and incoherent dynamics of the gas in the 
trap is for our toy model quite generally solved by 

(0n)(i) - (^„)(0)e-'(<-'^)*/''e-r"*/2 (3.108) 

and 

Nn{t)=N„{0)e-^-'+N„{l-e-^"*) , (3.109) 

respectively. In the limit t — > oo, the average of the annihilation operators 
{(j)n)(t) thus always vanishes but the average occupation numbers Nn{t) re- 
lax to the equilibrium distribution iVn = (e'^^'^""''^ — 1)~^. Although this 
appears to be an immediately obvious result, its importance stems from 
the fact that it is also true if we tune the potential energy bias AV°'^ such 
that at low temperatures the groundstate Xo(x) acquires a macroscopic 
occupation, i.e., A^o ^ 1. The gas therefore never shows a spontaneous 
breaking of the U{1) symmetry, in agreement with the notion that we are 
essentially dealing with an ideal Bose gas in the grand-canonical ensem- 
ble PQ] . The reason of the absense of the spontaneous symmetry breaking 
can also be understood from our stationary solution of the Fokker-Planck 



equation in equation (3.104), which shows that the probability distribution 
for |0o| is proportional to the Boltzmann factor e~^^'^o~t^)\'i'o\ [j^ i]^q ^jg- 
generate regime of interest. Hence, the corresponding Landau free energy 
Fl{\<Po\) = (eq ~ m)I'/'oP never shows an instability towards the formation 
of a nonzero average of |0o| due to the fact that eg — /i can never become 
less then zero. Once we introduce interactions between the atoms in the 
gas, this picture fully changes. 
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3.4 Condensate formation 

We are now in a position to discuss an interacting Bose gas, because in 
an interacting Bose gas, the gas is roughly speaking its own thermal reser- 
voir. The Fokker-Planck equation therefore turns out to be quite similar 



to the one presented in equation (3.103). To be more precise, however, we 



should mention that we aim in this section to describe the formation of a 
condensate in an interacting Bose gas under the conditions that have re- 
cently been realized in experiments with atomic rubidium, lithium, sodium 
and hydrogen gases. In these experiments the gas is cooled by means of 
evaporative cooling. Numerical solutions of the Boltzmann equation have 
shown that during evaporative cooling the energy distribution function is 
well described by an equilibrium distribution with time dependent temper- 
ature T{t) and chemical potential /i(t), that is truncated at high energies 
due to the evaporation of the highest energetic atoms from the trap [ pi[ . 
Moreover, in the experiments of interest the densities just above the critical 
temperature are essentially always such that the gas is in the weak-coupling 
limit, which implies in this context that the average interaction energy per 
atom is always much less than the energy splitting of the one-particle states 
in the harmonic trapping potential. 



3.4.1 Weak-coupling limit 

For a trapped atomic gas the appropriate one-particle states are the eigen- 
states Xn(x) of the external potential V'^^{x). As a result the Fokker-Planck 
equation for the interacting gas will now quite generally be given by 






P[0*, (/.; t] = -J2^\ ((en - Kt))Sn^n' + ^iS^) </-„' (3.110) 

n,n' I 

m,m' J 

n.n' ^ I 



m,m' 



^E7U^^^^^n'+ E K^n^n'^^'C-^m'^P^^;*] 



2 

n.n 
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For this equation to be useful, however, we need to determine the various 
selfenergies and interactions. To do so, let us for simplicity first consider 
the normal state of the gas which is usually in the weak-coupling limit 
|?iE^^, I ^ |en — fn'l for the experiments with atomic alkali gases that are 
of interest here. In this regime we can neglect the nondiagonal elements 
of the sel fenerg ies and also the nonlinear terms in the right-hand side of 
equation ( [3.110 ). The Fokker-Planck equation then reduces to 






dt 



(3.111) 



d 









with hTin — fi-Sn^n (e^ — fi) and e^ — en + Sn- This equation looks 
identical to the one we derived for our Caldeira-Leggett toy model but 
describes in fact quite different physics because now the selfenergies are not 
due to an interaction with a reservoir but due to the interactions between the 
atoms themselves. This is, however, completely hidden in the selfenergies as 
we show now explicitly. Performing the same manipulations as in sections 



2.5 and 2.8.2, we find in the many-body T-matrix approximation first of all 
that 



?iEW = (ns(,-)) * - 2 ^ vW;„,m(4 + C - ^^^)N„ 



(3.112) 



de 






' n.m;n,m 



W- 



^- - (e'n + e'm - 2m) 



and similarly that 



m 
+ 2l] V„1m;n,m(4 + ^^^ - 2^)(1 + 2N„,) 



(3.113) 



Our next task is therefore the evaluation of the various many-body T- 

n,m;n',m' 



matrix elements K!t!:^'^, (e) ^ {T^^^^,%, (e) + Ti+^l<;^, {e))/2 that occur 



in these expressions, which is easily carried out with our experience of the 
homogeneous gas and by noting that the Lippmann-Schwinger equation is 
in fact an operator equation that only needs to be expressed in a new basis. 
In terms of the appropriate matrix elements of the interatomic interaction. 
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I.e., 

^„,m;„'.m' ^ J d^ J d^' X.l{^)X*mi^')V {^ ~ ^')Xn' {^)Xn.' {^') , (3.114) 

the retarded and advanced components of the many-body T matrix thus 
obey the Bethe-Salpeter equation |9^ 

^n,m;n',m'(^) = Vn.min'.m' (3.115) 

, Y- . . 1 + Nn" + iVm" ^(±) , . 

+ 7. ''n,m;n" ,m" 4- // ., r, , x -^ n".m":n',m' l^J 

and the optical theorem 

Ti'^ln'.n.'i^) - TL:i,n',m'i^) = "2^* E '^(^ " (^n" + C' " 2^.)) (3.116) 

X 7^^:'J.;n",m"(f)(l +A^n" +iV^")r^"L';n',m'(^) • 

Finally, wc also need 

rn1m;„',m'(e) = -27rz ^ 5(6 - (6'„„ + e^^,, - 2^.)) (3.117) 

X -V,m;n",m"(^)-'^n"-'^m"^n",m";n',m'(^) 

and the Keldysh component 

T^,m;n',m'{^)^-27Tl ^ ^(^ " (4" + C" - 2^)) (3.118) 

n'^.m" 

X r^:'2;„",n."(e)(l + A^n" + N^„ + 2N^„ N^„)T^,]^„ .^^, ^^,{e) . 

Notice that the appearance of a (5 function in the last three equations in prac- 
tice never leads to any problems, because in the experiments with trapped 
atomic gases the number of particles is always so large that the critical 
temperature of the gas is much larger than the energy splitting of the trap 
and the discrete sum over states is well approximated by an integral over a 
continuous spectrum. 

Making use of the above relations we find that for the states with thermal 
energies the rcnormalized energy is given by the pseudopotential result 

4(^) - en + -^^^— / rfx x;(x)n(x, i)xn(x) , (3.119) 

m J 

with n(x, t) the density profile of the gas cloud. For the states near the 
one-particle ground state there is, however, a correction given by 

2E E i<^^-.^-""'^-'"^-'" .,(,)+cw4w-c-w 



m n",m" 
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and V;(+i^„..„..(i) = K!2;n':,m",(4"(0 + e'^4t) - 2^l{t)). In the weak- 
couphng limit, the renormahzation of the one-particle energies plays by 
definition not such an important role and we are much more interested in 
the imaginary part of the retarded and advanced selfenergies. They are 
given by 

i?„(t) = 27r^ Y, S{e'^„{t) + e'^„{t)~e'^{t)-e'^{t)) (3.120) 

m n",m" 

X lK!in",™"(i)l'[(l + N^"it) + N^„{t))N^{t) - iVn" W^m" W] • 

Furthermore, the Keldysh component of the selfenergies equals 

aE^(i) = -4vr^ Y. -5(4" W + C'(i) - 4W - e'mit)) (3.121) 



m n",m" 



As a result the quantum Boltzmann equation for the gas becomes 

l^n(t) - -TT{t)N„it) + C(t)(l + N^m , (3.122) 

with the rate to scatter out of the state Xn (x) given by 

4-7r 

rrw^Y^I E '5(e:,"W+C'W-4w-C(i)) (3.123) 

m n",m" 

x\V^!l;r."^m"it)\'{l+N^"imi+N^„it))N^it) 

and the rate to scatter into this state by 

rn(i) = YE E ^(4"(i)+e:n"(0-4(0-£:.(0) (3.124) 

m n",m" 

x|K!rf,m"(i)l'A^n"(t)iVm"(i)(l +A^m(i)) , 

as expected from Fermi's Golden Rule. 

The quantum Boltzmann equation found above, in combination with the 
expressions for the shift in the energy levels, fully describes the dynamics 
of the gas in the normal state. However, close to the critical temperature 
when the occupation numbers of the one-particle groundstate Xo(x) start 
to become large, we are no longer allowed to neglect the nonlinear terms 
in the Fokker-Planck equation in equation ( ^.110 ). Applying again the 



Hartree-Fock approximation, which is essentially exact in the weak-coupling 
limit as we have seen in the equilibrium case, and substituting the ansatz 
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P[(/)*,(/);i] = Po[(l)o,(l)o;t]Pi[(j)'*,(j>';t] into equation ( |3.110|) , we obtain for 
the probabiUty distribution of the condensate 

ih^^Po[ro,<t'o;t] (3.125) 

= -^ (eo + ?i4^^ - ^(t) + r^+o;0,oWI'^op) cj,oPo[<l>*o, '/'o;i] 

1 92 



PqOcPo 



-?iS^Po[0S,0o;i] , 



where it should be remembered that for consistency reasons the selfenergies 
need to be evaluated at the energies eg(i) — e'Q{t) + Tq o(^)l<^oP- In 
equilibrium the fluctuation-dissipation theorem then reads 

«^n = -^?iS^(e^-/i + T^So.ol0o|') , (3.126) 

which guarantees that the probability distribution for the condensate relaxes 
to the correct equilibrium distribution 



-fo[0S' 0o; 00] ex exp < -/3 Uq - /i 



-'0, 0:0,0 




The Fokker-Planck equation, and hence the Boltzmann equation, for the 
noncondensed part of the gas is the same as in the normal state because for a 
discrete system the contributions from the condenstate to the selfenergies is 
already included in our expressions for hJ^n and needs not be separated 
out explicitly. This implies of course, that if we evaluate the rates r°"*(i) 
and rj^' (t) by replacing the sum over states by an integral over a continuum, 
we first need to separate the contribution from the condensate to obtain the 
correct result. 

In summary we have thus found that in the weak-coupling limit the 
dynamics of Bose-Einstein condensation in an external trapping potential is 
described by the coupled equations in equation ( 3.122 ) and equation ( ^.125 ), 



i.e., by a Boltzmann equation for the noncondensed part of the gas and a 
Fokker-Planck equation for the condensate, respectively. If we apply these 
equations to a Bose gas with a positive scattering length, we essentially 
arrive at the theory recently put forward by Gardiner and coworkers [ p3[ . 
Indeed, neglecting the energy shifts and assuming that the noncondensed 
cloud is in equilibrium with a chemical potential /i > eo, our Fokker-Planck 
equation for a U{1) invariant probability distribution of the condensate 
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Fig. 17. Typical evolution of the condensate probability distribution during Bose- 
Einstein condensation. The slowest relaxation rate is 1/ro ~ 5.7. The inset shows 
the evolution of the average condensate number and its fluctuations. 



becomes 

/3„^ d 



.(+) 



eo-fi + n:o',o,o\M']<PoPo[\M;t] (3.127) 



^^0-. 



HE, 



K 



^0 

32 



.(+) 



eo-^l + UV..o.o\M')<f>oPo[\M■,t] 



-PQ[\M;t] 



2 u(PqU(Po 

It is equivalent to the Fokker-Planck equation for a single-mode laser, since 
introduction of the dimensionless time r — t(iE^/8)(2/3Tg q.q q)^/^^ t}jg 

dimensionless condensate number / = |0o|2(2/9Tg oo.o)^^^' ^'^d the so-called 

'pump parameter' a — 2/3(/i — eo)/(2/3Tg g q g)^/^ leads to the Fokker-Planck 

equation 

d 1 

' "'^ ^ (3.128) 



^^[^^^]-^P(^-«)^ + 4^a/ 



P I\t] 



"well-known from laser theory [ p4[ . The typical behaviour of P[/;t] during 
Bose-Einstein condensation in shown in figure 17. 
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Calculating (/)(t) from equation ( 3.12S| ) by making use of the mean- 



field assumption {P){t) — ((/)(t))^, precisely reproduces the theory of 
Gardiner et al. A comparison with the experiments of H.-J. Miesner et al. 
shows that the theory of Gardiner and coworkers is unfortunately not fully 
satisfactory l^^jQ- The reason for this is presumably that these experi- 
ments are really in the strong-coupling limit, where the above theory does 
not apply. As mentioned previously, the theory is however applicable to 
gases with a negative scattering length. Since we know from our experience 
with homogeneous gases that the chemical potential will never be larger 
than e'fj in this case, we clearly see from the equilibrium probability distri- 
bution Pq [(J3q , 00 ; oo] that a condensate in a gas with effectively attractive 
interactions is only metastable and will ultimately collapse to a dense state 
after it has overcome an energy barrier, either by thermal fluctuations or 
by quantum mechanical tunneling. This is clearly in full agreement with 



the discussion in section 3.1, Furthermore, we now also understand how 
the dynamics of the collapse must be treated. It requires the solution of 
a stochastic, dissipative nonlinear Schrodinger equation, in which the dis- 
sipative terms and the noise obey an appropriate fluctuation-dissipation 
theorem. Having arrived at this important result, we are finally left with a 
discussion of the strong-coupling limit. 

3.4.2 Strong-coupling limit 

The difficulty with the strong-coupling limit is that the appropriate one- 
particle states of the gas are no longer the states Xn(x) of the trapping 
potential because these states are now strongly coupled by the mean-field 
interaction. For our Caldeira-Leggett model this problem could be easily 
resolved by introducing the states Xn(^) that incorporate the average in- 
teraction with the reservoir. However, for a trapped Bose gas the same 
procedure would lead to an expansion into states Xn(^;^) that paramet- 
rically depend on time since the mean-field interaction depends on time 
also. In principle, we thus have to deal with the situation that not only the 
occupation numbers vary with time but also the states to which these oc- 
cupation numbers belong. This results in highly nontrivial dynamics since 
it is a priori clear that an adiabatic approximation, in which we neglect the 
nondiagonal matrix elements / dx Xni^jt){9/dt)Xn'{^'-,t)i is not valid, be- 
cause of the lack of a separation of time scales. Notice that we were able to 
circumvent the above problems for a Bose gas with effectively attractive in- 
teractions due to the fact that on the one hand the growth of the condensate 
essentially takes place in the weak-coupling regime and that on the other 
hand the collapse, which certainly occurs in the strong-coupling regime, re- 
quires in first instance only knowledge of the collisionless dynamics of the 
condensate. 

Notwithstanding the above remarks, we can make progress in the strong- 
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coupling limit if we realize that in this limit the typical length scale on 
which the (noncondensate) density of the gas varies is always much larger 
than the trap size. We are therefore justified in performing a gradient 
expansion. To do so in a systematic fashion, it is convenient to express the 
various selfenergies and interactions in the coordinate representation, i.e., 
by introducing 

?.I]W-^(x,x';6) = E^"W^^S;''W^n'(x') (3.129) 



and similarly 

T/W'^(x,x';y,y';6) (3.130) 

= E E x„(x)xm(y)K!S^,^,(^)x;'(x')x:„'(y') • 



n,n' m,m' 



In this coordinate representation, our Fokker-Planck equation acquires in 
all generality the form 

z?i^P[r , cj);t]^-J d^d^' j^ (3.131) 

{ ((-^ + ^°^(x) - A*(t)) <5(x- x') + as(+)(x,x')) 0(x') 
+ j dydy' F(+) (x, x'; y, y')0*(y)0(y')'^(x')| P[4>* , </-; t] 
dxax. 



/ dxdx.' 

2 



,50* (x) 
((-^ + V^-{^) - fi{t)^ 5(x- x') + ?iS(-)(x',x)^ r (x') 

+ ydydy'T/(-)(x',x;y',y)r(x')'^*(y')0(y)}^[r,0;i] 

W){'""(^'^'^ 

+ J dydy' y^(x,x';y,y')0*(y)0(y')|^r,0;t] • 

To study its physical content and to familiarize ourselves with how to per- 
form the desired gradient expansion, we first analyze again the normal state 
of the gas in which we can neglect the nonlinear terms representing the effect 
of the condensate. 

In the normal state the condensate is absent and we are only interested 
in the behavior of the quantity (0(x)(/)*(x'))(t). From the above Fokker- 
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Planck equation we see that it obeys the equation of motion 

z^^(0(x)<^*(x'))(t) (3.132) 

+ ;ii](+)(x,x"))(0(x")r(x'))W 
+ ;ii](-)(x",x')')Wx)<^*(x"))(i) 



-iftS^(x,x') . (3.133) 

Introducing now the Wigner distribution pJ7| ] by means of 

Af(x,k;i) + i= /dx' e-*-^'(0(x + x72)0*(x-x72))(i) (3.134) 
and similarly the selfenergies ?iE'^^)'^(x, k) by 

;iE(±)'^(x, k) = /" dx' e-*-^'?iE(±)'^(x + x72, x - x72) , (3.135) 

we find after substituting these relations in lowest nonvanishing order in the 
gradient 9/9x the expected quantum Boltzmann equation ||77|,|8i^] 

^iV(x,k;i) (3.136) 

, iMZElM . ±NU k-t) - iMZElM . Amx k-t) 
^h dk d^^^""'^'*' h ax dk^^""'^'*' 

= -r°"'(x, k; t)N{x, k; t) + r"(x, k; t)(l + iV(x, k; t)) , 

where the energies of the atoms have to be determined selfconsistently from 

e'(x,k;i) =e(k) + l^'='^(x)+Re[;iE(+)(x,k;e'(x,k;i)-^(t))] . (3.137) 

Moreover, the expressions for the retarded energy ?iS(+)(x, k; e'(x, k; i) — 
/i(i)), and the collision rates r°^*(x, k;i) and r'"(x, k;i) turn out to be 
identical to the ones for a homogeneous gas [£l|,|9|,^. We only need to 
replace the one-particle energies e'(k; t) and occupation numbers N{k; t) by 
their local counterparts e'(x, k;t) and A^(x, k;i), respectively. Notice also 
that the same substitution has to be made in the Bethe-Salpeter equation 
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for the many-body T matrix, which impUes that the effective interaction 
now impHcitly depends both on the position x and the time t. 

As an example, let us consider the gas not too close to the critical 
temperature. Then the effective interaction is well approximated by the 
two-body T matrix and we simply have 

.'(X, k; t) = e(k) + F-(x) + 5!I!!M^ . (3.138) 

m 

As a result the equilibrium solution A^(x, k; oo) of the quantum Boltz- 
mann equation equals a Bose-Einstein distribution evaluated at the en- 
ergy e(k) + V^°'^(x) -|- 87m(x)a?i /m — /x. The equilibrium density profile 
n(x) = J dk iV(x, k; cxD)/(27r)'^ thus corresponds to a local-density approxi- 
mation but is indeed very accurate for the atomic gases of interest. Apart 
from describing the relaxation of the gas towards equilibrium, the quantum 
Boltzmann equation can also be used to study the collective modes of the 
trapped gas above the critical temperature. Using the same methods that 
have recently been applied to inhomogeneous Bose gases by Zaremba, Grif- 



fin, and Nikuni | IOC ] , we can for example derive the hydrodynamic equations 



of motion the gas. Moreover, the coUisionless modes are obtained from 

f) Tile f) 

^Af(x,k;t) + — ~iV(x,k;i) (3.139) 

ot m ox 

Solving this equation is equivalent to the Generalized Random-Phase Ap- 
proximation |83| . More important is that the above coUisionless Boltzmann 



equation can be shown to satisfy the Kohn theorem [101| by substituting 
the shifted equilibrium profile iV(x — Xcm(i), k — {m / h)dKcm{t) / dt; oo). For 
an harmonic external potential the center-of-mass of the gas cloud then pre- 
cisely follows Newton's law d^x.cm{t) / dt^ — — VF'^''(xcm(0)7 which shows 
that the gas has three coUisionless modes with frequencies that are equal to 
the three frequencies of the external trapping potential. 



The full quantum Boltzmann equation in equation ( 3.136| ) also describes 



the relaxation of the gas towards equilibrium, ft can thus be used to study 
forced evaporative cooling by adding an appropriate loss term to account for 
the escaping high-energetic atoms. In agreement with our previous picture 
we then expect that, if the evaporative cooling is performed slowly, the 
distribution function iV(x, k;i) can be well approximated by a truncated 
equilibrium distribution with a time-dependent temperature and chemical 
potential. As a result, we can use the results from the homogeneous gas 
to show that in the critical regime the gas again becomes unstable towards 
Bose-Einstein condensation. This is the case when e'(x, 0;t) < /i(i) for 
positions near the center of the trap where the density of the gas is highest. 
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If that happens we are no longer allowed to neglect the nonlinear terms in 
the Fokker-Planck equation and some corrections are neccesary. Applying 
again the Hartee-Fock approximation, which is correct sufficiently close to 
the critical temperature where ksT ^ 47rno(x, i)a?i^/m, we find ffist of 
all that the quantum Boltzmann equation for the noncondensed part of 
the gas has the same form as in equation ( |3.136| ). We only have to add 
the by now familiar mean-field interaction 2no(x, i)T(+)(O,0) due to the 
condensate to the renormalized energy e'(x, k;i). Finally, the scattering 
rates r°"'(x, k; t) and r'"(x, k; t) have to be adjusted in precisely the same 
way as in the homogeneous case to account for the macroscopic occupation 



of the condensate |98 , 102, 103| 



In addition, the resulting quantum Boltzmann equation for the noncon- 
densed atoms is coupled to a Fokker-Planck equation for the condensate. 
After a gradient expansion the latter equals 

tn^Po[<f>*o,^o;t] (3.140) 



dt 



fd^ -^ (-%^ + y-(x) - M(i) + ?.sW(x,o) 

J 000 (x) V 2m 

+ r(+)(O,O)|0o(x)p') </)o(x)Fo['/'o,'/'o;i] 



5 ( fi^V^ 



''-m^A-^^^'''^^-^-^^'^^^'''''^-^'^ 



+ r(+)(o,o)|0o(x)n 0S(x)Po[0S,0o;t] 

4/-^- .0o(x)l5(x) ^^"^^'"^^°[^°'^°^^J' 

where ?il]'^-''^(x, 0) denotes ?iS(*^'^(x, k = 0). In more detail we have 
for the real part of the retarded and advanced selfenergies that S'(x, 0) = 
Re[^I](=^'(x,0;e"(x,0;t) - /i(i))] and e"(x,0;t) the local energy of a con- 
densate atom. Moreover, using our picture of the evaporative cooling of 
the gas, the imaginary parts follow most conveniently from the fluctuation- 
dissipation theorem 

zi?(x, 0) = -^ai]^(x, 0) (3.141) 

X f-^ + e'(x, 0; t) - fi{t) + r(+)(o, 0)\M^)\A 

\ 2m J 

where ?il]^(x, 0) = /il]'^='=^'^(x, 0; e"(x, 0; t)—fj,(t)) can in a good approxima- 
tion be considered as independent of |(/)o(x)|. Near equilibrium the Fokker- 
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Planck equation thus simplifies to 

ih-^Po[\M,t] (3.142) 



dt 



4 J ^<50o(x) V 2m 

+ T(+)(O,O)|0o(x)P') 0o(x)Po[|0o|;i] 

f /dxai]^(x,0)— ^f-^+e'(x,0)-M 

+ r(+)(O,O)|0o(x)p')0*(x)Po[|0o|;i] 
l/dx;.E^(x,0) ^^^^^g^^^^^ Po[|^o|;^] 



and the probability distribution of the condensate relaxes to the expected 
stationary state 

Po[|(/.o|;oo]cxexp|-/3y dx(/.*(x)f-^— +e'(x,0)-^ (3.143) 

+^:^i;iM)|,o(x)pVo(x) 



From this result we conclude that Eqs. ( 3.136| ) and (3.140) give an ac- 



curate discription of the dynamics of Bose-Einstein condensation in the 
strong-coupling regime. Unfortunately, the solution of these coupled equa- 
tions is not possible analytically and we therefore have to refer to future 
work for a numerical study and a comparison of this theory with experi- 
ment. We here only remark that the semiclassical approximation in this 
case amounts to the Thomas-Fermi result for the condensate density 

'^o(x; t) = 7^(^^J^^ (pit) - e'(x, 0; t)) e{^i{t) ~ e'(x, 0; t)) . (3.144) 

Such an approximation implicitly assumes, however, that the dominant 
mechanism for the formation of the condensate is by means of a coher- 
ent evolution and not by incoherent collisions. This might be expected 
not to be very accurate on the basis of the homogeneous result that the 
second mechanism only diverges logarithmically with the system size [ po[ , 



whereas the first diverges linea rly |104 |. N evertheless it seems that only a 



full solution of equations ( 3.136| ) and ( ^.14C| ) can settle this important issue. 



Les Houches summer schools 



3.5 Collective modes 

Although we have up to now primarily been focussed on the dynamics of 
Bose-Einstein condensation itself, it is interesting that we have in fact, sim- 
ilarly as for temperatures above the critical temperature, also arrived at an 
accurate discription of the collective modes of a condensed gas at nonzero 
temperatures. Since these collective modes are receiving much attention 
at present, we now briefly want to discuss how such a description emerges 
from our Fokker-Planck equation. Since most experiments at present are 
performed in the coUisionless limit, we consider only this limit here. It is 
in principle possible, however, to derive the hydrodynamic equations of mo- 



tions as well [100|. Let us for simplicitly first assume that the temperature 
of the gas is sufficiently far below the critical temperature that the many- 
body T matrix is well approximated by the two-body T matrix, but not so 
low that the Hartree-Fock approximation is no longer valid. In this regime 
we have that TiS^^'^^ (x, k) — STrn'iXjfjah^/m and the dynamics of the non- 
condensed part of the gas is again determined by the coUisionless quantum 
Boltzmann equation given in equation (|3.139), where the total density now 



consists of two contributions 

/dk 
---3 7V(x,k;i) . (3.145) 

The dynamics of the condensate follows again from the nonlinear Schroding- 
er equation 

*?i^((/)(x))(t) (3.146) 

^ + l^-(x) - A* + l!^(2n'(x, t) + no(x, t))\ (0(x))(i) , 

which can be rewritten as the continuity equation 

d 

— no(x, t) ^ - V • (v,(x, t)no(x, t)) (3.147) 

together with the Josephson relation for the superfluid velocity 

|^v,(x,t) (3.148) 

ot 

= -V fv^"'^(x) -/.+ ^^^(2n'(x,i) +«o(x,t)) + iTOv2(x,f) 
V ?Ti 2 



Here, we used that ((?!)(x))(t) = ^no(x,i)e*''('''*\ Vs(x,i) = hVe{x,t)/m, 
and have applied the Thomas-Fermi approximation, which is appropriate 
in the strong-coupling limit. It will presumably not come as a surprise that 
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these coupled equations for iV(x, k;t), no(x, i), and Vs(x, i) can be shown 
to fulfill the Kohn theorem exactly. This is an important constraint on any 
theory for the collective modes of a trapped gas, as stressed recently by 



Zaremba, GrifRn, and Nikuni |10C] 



As a first step towards this goal we here consider the time-dependent 
Hartree-Fock theory for the collective modes of the Bose condensed gas, in 
which the coUisionless Boltzmann equation in equation ([3.13S) is coupled 



to the nonlinear Schrodinger equation in equation (3.146|). It turns out 



that a numerical determination of the eigenmodes of the coUisionless Boltz- 
mann equation is rather difficult. Therefore, we have recently proposed a 



variational approach to solve these coupled equations |105]. Denoting the 
anisotropic harmonic trapping potential used in the experiments p6|,p7| by 
V'^^{x) = ^ Tnuj'^x'^/2, we first of all take for the condensate wave function 
the gaussian ansatz 

(0(x))(t) (3.149) 

^j (^ -rnqjit) dqj{t)\\ 




2q^{t) V fi dt J I ' 

with the functions qj(t) describing the 'breathing' of the condensate in three 
independent directions. Clearly, this is the anisotropic generalization of 



the gaussian ansatz that we used in section 3T to study the stability of a 
condensate with negative scattering length. Although it is not immediately 
obvious that such an ansatz is also appropriate in the strong-coupling regime 
of interest here, it has been shown by Perez-Garcia et al. [^ that at zero 
temperature the assumption of a gaussian density profile nevertheless leads 
to rather accurate results for the frequencies of the collective modes. We 
therefore anticipate that at nonzero temperatures the same will be the case. 
Next we also need an ansatz for the distribution function A'^(x, k;t) of 
the thermal cloud. Here, we consider the Bose distribution 



(3.150) 




ma;2 x] 



a 



it) 



where aj{t) are again three scaling parameters, N' — N — Nq is the number 
of noncondensed particles, (aj) denotes the equilibrium value of aj{t) that 
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Fig. 18. CoUsionless modes in a Bose condensed *7Rb gas. Tlie full lines give 
the results with a Bose distribution ansatz for the distribution function of the 
thermal cloud. Curves 1 and 2 correspond to the m = in and out-of-phase 
modes, respectively. Similarly, curves 3 and 4 give the frequencies of the m = 2 
in and out-of-phase modes. The experimental data is also shown. Triangles are 
for a m = mode and circles for a. m = 2 mode. 



in general is greater than 1 due to the repulsive interactions between the 
atoms, and C(3) — 1.202. The functional form of the distribution function 
is physically motivated by the following reasons. First, by integrating over 
the momentum ?ik we see that the time-dependent density profile of the 
noncondensed cloud n'(x, i) is a time- independent function of Xj/aj(i), as 
desired. Second, for such a density profile we must have that the local 
velocity h{kj){x.,t)/m is given by {xj / aj{t))daj (t) / dt. Third, we want the 
equilibrium distribution to be only anisotropic in coordinate space, and not 
in momentum space. This explains the appearance of the factor l/{aj)'^ in 
the exponent. 

Using the ideal gas result for the number of condensate particles and sub- 
stituting the above ansatz for the condensate wave function and the Wigner 
distribution function into the nonlinear Schrodinger and collisionlcss Boltz- 
mann equations, it is straightforward but somewhat tedious to derive the 
six coupled equations of motion for the 'breathing' parameters qj{t) and 
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aj{t). Linearizing these equations around the equihbrium solutions (qj) 
and (ttj), respectively, we can then also obtain the desired eigenmodes and 
eigenfrequencies. The results of this calculation for the trap parameters of 
Jin et al. |g^ are shown in figure 18 together with the experimental data. 
Note that the trap has a disk geometry with ujx = ^y ^ <^r and ujz = yS^r- 
Due to the axial symmetry, the modes of the gas can be labeled by the 
usual quantum number m that describes the angular dependence around 
the z-axis. On the whole there appears to be a quite reasonable agreement 
between this simple theory and experiment. The only exceptions are the 
two data point almost halfway in between the in and out-of-phase m = 
modes. A possible explanation for this discrepancy is that experimentally 



both modes are excited simultaneously |10(;]. To make sure of this, how- 
ever, we need to determine the oscillator strength of the various modes 
and it appears that the in-phase mode indeed dominates the response for 



temperatures larger than about 0.6Tc [105|. 

4 Outlook 

With the latter remark we end this course on the use of field-theoretical 
methods for the study of the equilibrium and nonequilibriuni properties of 
trapped atomic gases. Although we have in principle illustrated all the 
necessary tools for an ab initio treatment of these new quantum systems, 
many interesting topics still need to be considered in detail. We already 
mentioned such topics as Fermi-Bose mixtures, the dynamics of condensate 
collapse, the strong-coupling theory for the formation of a condensate in 
either atomic Bose of Fermi gases, and the hydrodynamics of single or mul- 
ticomponent atomic gases. Several other important areas of research, which 
can also be easily addressed within the context of quantum field theory, are 
atom lasers, the damping of collective modes, the behaviour of spinor con- 
densates, the dynamics of topological excitations such as kinks, vortices and 
skyrmions, quantum critical phenomena, the optical properties of superfluid 
gases, and two-dimensional phase transitions. In view of these and many 
other possibilities that are now experimentally feasible, it appears certain 
that the physics of degenerate atomic gases will remain very exciting for 
years to come. 
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